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Low-Cycle Fatigue of Reactor 
Structural Materials 


By C. E. Jaske, J. S. Perrin, and H. Mindlin* 


Abstract: Low-cycle fatigue and creep—fatigue are important 
considerations in evaluation and development of design in- 
formation for reactor structural materials, Strain-controlled 
low-cycle fatigue tests of axially loaded specimens of Alloy 
800, 304 stainless steel, and 316 stainless steel have been 
conducted at several laboratories, Fatigue curves for these 
three alloys have been developed at temperatures from 70 to 
1500°F. Hold times at peak strain and slower strain rates have 
been shown to greatly reduce cyclic-fatigue resistance at 
temperatures between 1000 and 1400°F. The linear-life- 
fraction damage rule has been used to evaluate creep—fatigue 
damage interaction in tension hold-time tests, but it does not 
provide a proper explanation of the damage mechanisms, These 
materials exhibited significant cyclic strain hardening at all 
temperatures studied, The hardening mechanism was in- 
fluenced by strain rate, number of strain cycles, and prior heat 
treatment of the alloy, Aging of material at test temperature 
prior to test had a significant influence on both creep—fatigue 
resistance and cyclic hardening at temperatures near 1100 to 
1200°F but had little effect at temperatures near 1000 to 
1050°F. 


During normal operation of a nuclear power plant, the 
pressure vessel and such related components as piping 
are subjected to transient loading during startup, 
shutdown, or variations in the steady-state operation 
conditions. During these transient conditions, com- 
ponents at elevated temperatures are subjected to 
nonuniform thermal stresses large enough to exceed 
the elastic limit and therefore to cause local plastic 
flow and introduction of residual stresses. The result of 
this can be damage from both fatigue and creep such 





*The authors are Senior Researcher, Associate Fellow, and 
Division Chief, respectively, at Battelle Memorial Institute, 
Columbus Laboratories, Columbus, Ohio 43201. 


that the component may fail in a relatively low number 
of cycles. 

Low-cycle fatigue is normally associated with the 
presence of cyclic macroscopic plastic strains. For 
many structural materials that are cyclically loaded 
under such conditions, fatigue failures occur in rela- 
tively short lifetimes of less than 10° cycles. The 
problem of low-cycle fatigue failure has received 
increasing consideration in recent years. The ASME 
Boiler and Pressure Vessel Code, Section III,' of the 
American Society of Mechanical Engineers (ASME) 
requires a fatigue analysis of pressure-vessel and related 
components used in nuclear reactors and provides 
fatigue design curves for use at temperatures up to 
700°F for carbon, low-alloy, and high-tensile steels and 
at temperatures up to 800°F for austenitic steels, 
Ni—Cr—Fe alloy, Ni—Fe—Cr alloy, and Ni—Cu alloy. 
Design rules for use of 304 and 316 stainless steels at 
temperatures above 800°F are covered? in ASME Code 
Case 1331-5. Rules for use of solution-annealed Alloy 
800 are to be added to the next revision of this Code 
Case. 

Before design rules for the above-mentioned 
materials can be developed, experimental information 
on their low-cycle fatigue behavior must be obtained. 
This review will cover the special experimental pro- 
cedures normally used in low-cycle fatigue testing. In 
addition, the types of low-cycle fatigue results which 
have been and are currently being obtained by investi- 
gators on three specific materials (Alloy 800, 304 
stainless steel, and 316 stainless steel) will be presented 
and discussed. These materials are all of current 
interest for nuclear reactor applications. 


REACTOR TECHNOLOGY, Vol. 15, No. 3, Fall 1972 








186 LOW-CYCLE FATIGUE OF MATERIALS 


EXPERIMENTAL PROCEDURES 


Low-cycle fatigue testing techniques have evolved 
relatively recently; a review that contains descriptions 
of approaches used at a number of laboratories has 
been published.* During low-cycle testing, specimens 
are normally run under strain control rather than load 
control. Usually specimens are tested using uniaxial 
loading, although some work is done by loading the 
specimen in bending. 

The most versatile type of machine for low-cycle 
fatigue testing is the closed-loop, electrohydraulic, 
servo-controlled testing system.* With such a system, 
any type of strain—time waveform desired can be 
programmed. For continuous-cycling tests, strain is 
usually programmed to follow a fully reversed tri- 
angular waveform (Fig. 1) with a constant strain rate, 
€, and a constant total strain range, Ae,. Initial loading 
is usually compressive (as shown in Fig. 1) because it is 
more representative of actual service conditions. Time- 











Oesias t 
g cS Cs / 
S ios aot Time 


Fig, 1 Schematic of program waveform for fully reversed 
continuous-cy cling fatigue test. 


dependent or creep deformation is usually introduced 
through a hold time, fj, at peak strain as illustrated in 
Fig. 2. At elevated temperature and significantly high 
strain levels, stress relaxation usually takes place during 
the hold period. Load is applied from a hydraulic 
actuator to the test specimen, and load is measured 
using a load cell. 

Two common types of specimen configuration 
(Fig. 3) are used in uniaxial loading tests; both have a 
circular cross section. The cylindrical type [Fig. 3(a)]| 
has a constant cross section and a longitudinal gauge 
section, whereas the hourglass type [Fig. 3(b)] has a 
minimum diameter at the center of the test section and 
a diametral gauge length. Other types include hollow, 
tubular specimens and specimens with a notched gauge 
section to localize failure. 
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Fig, 2 Schematic of program waveforms for fatigue tests with 
hold times. 
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(b) Hourglass Shaped 


Fig. 3 Typical configurations of low-cycle fatigue specimens, 


Strain can be controlled and measured with either a 
longitudinal (Fig. 4) or diametral (Fig. 5) extensome- 
ter. The longitudinal extensometer is used with a 
specimen of uniform gauge section and has the 
advantage of measuring the axial strain directly. How- 
ever, the extensometer must be mounted carefully so 
that the specimen surface is not disturbed at the point 
of attachment. Also, the specimen must be machined 
with very small diameter variations along the gauge 
length, and temperature must be uniform along the 
gauge length. The diametral extensometer has the 
advantages that temperature has to be controlled over a 
short length, the extensometer is following the be- 
havior of the specimen right at the location of ultimate 
failure, and the possibility of buckling failure is 
minimized. However, to control axial strain with an 
hourglass-shaped specimen requires that the diametral 
strain and load be automatically combined and con- 
verted to equivalent axial strain during testing.° 

Heating the specimen to the desired test tempera- 
ture can be accomplished by such methods as direct 
resistance heating of the specimen, induction heating, 
or putting a furnace around the entire specimen. 
Temperature is normally measured by attaching 
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Fig.4 Apparatus for low-cycle fatigue testing of cylindrical 
specimens at elevated temperature. 





Fig.5 Apparatus for low-cycle fatigue testing of hourglass- 
shaped specimens at elevated temperature. 
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thermocouples directly to the specimen surface. In 
some cases infrared pyrometry is also used. 


MATERIALS 


In recent years a number of experimental programs 
have been conducted under the sponsorship of the 
U.S. Atomic Energy Commission (USAEC) to evaluate 
the low-cycle fatigue resistance of structural reactor 
materials. An extensive investigation of the plastic 
fatigue properties of 304, 348, and 316 stainless steels 
was conducted at Nuclear Systems Programs, General 
Electric Company (GE-—NSP), from 1964 to 1969 
(Ref. 6). Constant-amplitude, axial-load,  strain- 
controlled tests of all three alloys were conducted’ at 
800, 1200, and 1500°F with strain rates of 4 x 10°, 
4x 107%, and 4x 10° sec!. The influence of hold 
time at peak strain was studied for 304 and 316 
stainless steels at 1200°F (Ref. 8). The notched low- 
cycle fatigue behavior of carbon steel, 2/, Cr—1 Mo 
steel, and 304 stainless steel was investigated at room 
temperature and 550°F in a 4-year program conducted 
at the Materials and Processes Laboratory of the 
General Electric Company.’ In another program con- 
ducted at GE—NSP, the constant-amplitude, axial-load, 
strain-controlled low-cycle fatigue behavior of Alloy 
800 was examined at temperatures!® from 800 to 
1400°F. In subsequent work at Battelle’s Columbus 
Laboratories (BCL), low-cycle fatigue studies were 
conducted on specimens from four different heats of 
solution-annealed Alloy 800 (Ref. 11). Best-fit fatigue 
curves were developed at 70, 800, 1000, 1200, and 
1400°F (Ref. 11), cyclic-stress—strain data were ob- 
tained at the same temperatures,!? and the influence 
of hold times at peak strain was examined at 1000, 
1200, and 1400°F (Ref. 13). 

In current work at BCL,'*"7 low-cycle fatigue 
design information on 304 and 316 stainless steels is 
being developed for use at temperatures up to 1200°F. 
This work includes determination of monotonic- and 
cyclic-stress—strain response, continuous-cycling 
fatigue resistance, the effect of aging and hold time on 
fatigue life, cumulative damage under variable- 
amplitude loading, and creep—fatigue damage interac- 
tion. Developing an understanding of the microstruc- 
tural phenomena involved in low-cycle fatigue failure 
of 304 and 316 stainless steels in the 800 to 1200°F 
temperature regime is the goal of experimental work 
being performed at Argonne National Laboratory 
(ANL).'*?7 The type of failure mode has been 
qualitatively related to creep—fatigue damage interac- 
tion, and the importance of material condition (i.e., 
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annealed, stress relieved, or aged) prior to testing has 
been emphasized. Also, fatigue crack-propagation ex- 
periments and computerized simulation of cyclic- 
loading-behavior studies are being initiated for both 
materials. Primary effort of a program at Aerojet 
Nuclear Company (ANC) has been directed to the 
influence of irradiation on creep—fatigue of 304 and 
316 stainless steels and Alloy 800 (Refs. 28-31). 
Irradiated and unirradiated specimens of both stainless 
steels have been tested at 1100°F under constant- 
amplitude, strain-controlled axial loading with and 
without hold times at peak tensile strain. The effect of 
irradiation on fatigue and creep—fatigue of these two 
stainless steels and of Alloy 800 at 1292°F was also 
examined. 

Details of the work summarized above are reported 
in the following sections for each of the three materials 
of primary interest: (1) Alloy 800, (2) 304 stainless 
steel, and (3) 316 stainless steel. 


Alloy 800 


Curves of total strain range vs. fatigue life are 
available at a strain rate of 4x 10% sec’ for mill- 
annealed Alloy 800 at 800, 1000, and 1100°F and for 
solution-annealed Alloy 800 at 1000, 1200, 1300, and 
1400°F (Ref.10). The fatigue resistance of this 
material decreased with increasing temperature, and 
the fatigue resistance of the mill-annealed material was 
superior to that of the solution-annealed material at 
temperatures below 1100°F. A few specimens tested at 
a strain rate of 4 x 10% sec” had slightly lower fatigue 
strength than comparable ones tested at the faster rate. 

Since the solution-annealed material was of pri- 
mary interest for use at high temperatures, subsequent 
studies were conducted on the same heat of material 
used by Conway! ° (specimens from this heat were 
designated as Group A material) and from three addi- 
tional heats of Alloy 800 (specimens were designated 
as belonging to Groups B, C, and D). Groups A and D 
were from standard heats of Alloy 800, whereas Group 
B was duplex solution-annealed and Group C had a 
high carbon content (0.09% C compared to 
0.05—0.06% C for the other heats). Best-fit fatigue 
curves for all the results on solution-annealed material 
at 4x 10° sec’! are shown in Fig. 6. Values of plastic 
strain range, Ae,, and elastic strain range, Aeg, were 
related to fatigue life, Ny, by the following polynomial 
equation: 


log Ae, (or log Aee) = 
A, +A, logN¢+A3 (logNy = (1) 
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where, A,, A, and A; are best-fit regression coeffi- 
cients. (See Table 1.) The value of Ae; (the total strain 
range) is simply the sum of Ae, and Aée for any 
particular value of Ny. For all cases but Ae, at 800°F, 
the data fell close to a straight line on the log-log plots 
and the constant A3 was deleted. Results for all four 
heats of material were about the same, and the fatigue 
results at each temperature were well approximated by 
a single Ae, vs. Ny curve. 
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Brinkman et al.?! have shown that irradiation of 
Alloy 800 to fluences of 0.3 to 5.0 x 107? neutrons/ 
cm? (E>0.1 MeV) at 1292 to 1382°F produced a 
significant reduction in fatigue life at 1292°F. For Ae; 
from 0.003 to 0.024, fatigue life was reduced by 
factors as large as 35. (See Fig. 7.) This severe 
degradation in fatigue resistance was caused by irradia- 
tion-induced embrittlement of the Alloy 800 and was 
correlated with corresponding changes in tensile 
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Fig. 6 Fatigue curves for solution-annealed Alloy 800 at a strain rate of 4x 10°? sec’. (From 


Conway’° and Jaske et al,’") 
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Table 1 Summary of Best-Fit Regression Coefficients for 
Solution-Annealed Alloy 800 Fatigue Curves 





Best-fit value for 


a Refi : 
indicated constant eferences from 
































Temp., Type of which data 
oF strain range AL, As A, were taken 
70 Aép 0.225 -0.499 11 
ee -1.76 0.177 
800 Aen -1,54 0.294 —0.124 11 
Mee 2.06 ~0.120 
1000 Aep 0.0143 —0.621 10, 11 
Mee -1.93 ~0.128 
1200 Aép 0.0421 —0.657 10; 11 
Aée 1.89 —0.0148 
1400 hep 0.264 —0.793 10, 11 
Mee 2.11 ~0.127 
10° Pie T T wees) summarized as shown in Fig. 8. At each temperature 
, Anneale ‘ ; 
0, Thermal contro! (B) 9240 hr at 1382°F the Ae, vs. Ny relation can be approximated by a 
0, Thermal control (A) 1541 hr at 1382°F ae soos, et ta Oe a > . P E 
- ©, irradiated (A) 3-9 x102" neutrons/cm2 (E >0.1 MeV) straight line for a given length of hold time , and these 
4 B, Irradiated (B) 1~5 x 1022 neutrons/cm? (E >0.1 MeV) straight lines all intercept the one-cycle-to-failure line 
Ww yo Ce ee re at the point which is a linear extrapolation of the 
oO —, Universal-slopes equation, average life 0.2 N¢ ; bee 
4 (),Fluence x 10722 no-hold-time fatigue curve. These results imply that 
r . . . . 
2 tensile hold times have a more detrimental influence o1 
© _ Control (B) fatigue as their length increases and that the detri- 
eres a mental effect is more pronounced at low strain ranges 
“10 I 8 
< es than at high strain ranges. Values of the slopes and 
e intercepts of the curves in Fig. 8 are summarized in 
Table 2. 
Irradiated (B) : , - 
13 zal a ae eal ; The hold-time results discussed above were all for 
10! 107 10° 104 10° 


FATIGUE LIFE (N¢), cycles 


Fig. 7 Effect of irradiation on the fatigue life of Alloy 800 at 
1292°F and at a strain rate of 8 x 10% sec’. (From 
Brinkman et al.°') 


properties using the modified method of universal 
slopes.2?°3? As shown in Fig. 7, the modified uni- 
versal-slopes equation provided a good prediction of 
the fatigue life of the unirradiated control specimens 
and the specimens irradiated to low fluences. For 
specimens irradiated to high fluences, it gave a con- 
servative prediction. 

Tension hold times [Fig. 2(a)] detrimentally affect 
the fatigue life of solution-annealed Alloy 800, whereas 
compression [Fig.2(b)] or combined tension and 
compression holds [Fig. 2(c)] have only a small effect 
or no significant effect on fatigue life.'* Results of 
tests at constant lengths of tensile hold time can be 
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material from Group A. Results of a few tension 
hold-time tests conducted on specimens from Groups B 
and C showed more heat-to-heat variation than for 
continuous cycling. Duplex solution-annealed material 
(Group B) had slightly inferior creep—fatigue resistance 
to the standard material at 1000 and 1200°F and 
about the same at 1400°F. High-carbon-content 
material (Group C) had markedly better creep—fatigue 
resistance than the standard material at 1000 and 
1200°F but was slightly worse at 1400°F. 

In recent work at BCL,'” 10-min tension hold-time 
tests were conducted on both as-received and aged 
(150 hr at 1200°F and 1000 hr at 1200°F) specimens 
from the Group D material. Tests were conducted at 
total strain ranges of 0.02 and 0.01 at 1200°F, and 
results were compared with previous data on Group A 
material.’ As shown in Fig. 8(b) (compare open, 
partially filled, and completely filled diamond 
symbols), greater hold-time fatigue resistance - was 
observed for the aged material than for the as-received 
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Fig. 8 Effect of hold time on fatigue life of solution-annealed Alloy 800 at a strain rate of 4 x 10°° 
sec! , (From Conway! ° and Jaske et al.’ *) 


Table 2 Slope and Intercept Values That Define Tension 
Hold-Time Fatigue Curves for Solution-Annealed Alloy 800 





Length of tension 


Slope of Ac; References from 








Temp., Intercept at hold time, vs. which data 
“F Ne= min Nf curve were taken 
1000 0.48 0 -0.47 10 
10 —0.50 13 
60 —0.57 13 
300 —0.71 13 
1200 0.57 0 —0.54 10 
10 —0.68 13 
60 —0.79 13 
300 —0.83 13 
1400 0.84 0 —0.63 10 
10 —0.74 13 
60 —0.83 13 
300 —0.88 13 
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material. As-received specimens had about the same 
hold-time fatigue lives as expected for the previously 
studied heat. However, specimens that were aged 
150 hr had about 14, times as much fatigue life, and 
those which were aged {000 hr had over 21, times as 
much life. It was also noted that increased hold-time 
fatigue life was accompanied by decreased cyclically 
stable stress range. Thus the aged material was sub- 
jected to less time-dependent damage than the as- 
received material because it had lower stabilized stress 
values. 

Examination of the cyclically stable relaxation 
curves from the tension hold-time tests showed that 
they could be represented by using the following 
relation of Gittus:*4 





In (00/0) = —* l+m (2) 


+m! 


where 69 = initial stress, ksi, at the beginning of the 
hold period 
o = instantaneous stress, ksi, at some point 
during the hold period 
t = time, min 
A =a constant 
m =a constant 


Using Eq. 2 to define stable stress-relaxation 
curves, creep—fatigue damage interaction was com- 
puted using a linear damage rule of the form suggested 
by ASME Code Case 1331-5 (Ref. 2). Total damage, D, 
was defined as 


D=iXn/N¢ + LAt/t, (3) 


where 1 = number of cycles applied at a particular load 
condition 
Ny = number of cycles to failure for continuous 
cycling at the same load condition 
A; =increment of time that creep stress is at a 
given level 
ty =stress-rupture life, hr, at the above- 
mentioned level of creep stress 


Results of creep—fatigue damage calculations that 
were made for Alloy 800 are presented in Fig. 9. 
Considerable scatter in the results is related to large 
variations in the values of the creep damage term 
(ZAt/t,) because ft, is very sensitive to the level of 
stress at these temperatures. Thus small errors in 
measurement of stress and normal scatter in rupture 
life will inherently cause large variations in creep 
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Fig.9 Creep—fatigue damage interaction for tension hold- 
time tests of solution-annealed Alloy 800. (From Jaske 
et al,!') 


damage summations. However, the general overall 
trend of the damage results was significant. Data were 
reasonably close to the D=1 curve when fatigue 
damage (2n/N¢) was greater than 0.5. At lower values 
of fatigue damage than 0.5, the results were signifi- 
cantly greater than D = 1. The damage values tended to 
increase as the length of tension hold time increased as 
shown in Fig. 10. 
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Fig. 10 Effect of length of tension hold time on total damage 
for solution-annealed Alloy 800. (From Jaske et al.'') 
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Incremental step tests have been used to obtain 
comparative monotonic and _ cyclic-stress—strain 
curves'? for solution-annealed Alloy 800. (See 
Fig. 11.) Results from the step tests provide a fairly 
good correlation with those from constant-amplitude 
tests. (Compare solid curves with plotted symbols in 
Fig. 11.) At strain amplitudes lower than the maximum 
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used in the step test, there is some deviation in results 
between the two methods of defining stable stress— 
strain behavior. This deviation indicates that prior 
history has an influence on such behavior; thus, in the 
variable-amplitude step test, the alloy has a “memory” 
for high-level strain cycles that influence subsequent 
strain cycles at low strain levels. Since the cyclic- 


Constant Amplitude, Continuous Cycling 
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Fig. 11 Comparison of monotonic- and cyclic-stress—strain curves for solution-annealed Alloy 800 at 


a strain rate of 4 x 10°° sec! (From Jaske et al,''’'?) 
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stress—strain curves fell significantly above the 
monotonic ones, a large amount of cyclic hardening 
occurred. Comparing the cyclic with the monotonic 
yield strength (Table 3) shows that the alloy cyclically 
hardened by factors from 2.3 to 4.4 with a maximum 
amount of hardening at 1200°F. 

The cyclic-stress—strain curves in Fig.11 were 
approximated by the following function: 


Ao/2 = K'(Aep/2)" (4) 


where Ao/2 = cyclically stable stress amplitude, ksi 
Aép/2 = cyclically stable plastic-strain amplitude 
K' = cyclic-strength coefficient, ksi 
n' = cyclic-strain-hardening exponent 


Knowing the modulus of elasticity, F, values of Ae, /2 
were defined by the relation 


Mey /2 = Aez/2 — Ao/2E (5) 


Analogous equations were used for the monotonic- 
stress—strain curves, and both monotonic- and cyclic- 
strength coefficients and strain-hardening exponents 
are summarized in Table 3. 

Data from constant-amplitude tests?! of Alloy 800 
showed that irradiation had no significant effect on 
cyclic-stress—strain behavior at 1292°F. (See Fig. 12.) 
The small increase in yield strength (20 to 23 ksi) from 
the monotonic to the cyclic case indicated only a small 
amount of cyclic hardening for this heat of material. 
The main differences between this material and that 
used in the other studies! ? are that the carbon content 
was 0.03% instead of 0.05 to 0.06% and that the final 
anneal was at 1904°F instead of at 2100°F. Further 
studies are necessary to determine why this heat of 
material has relatively stable cyclic-stress—strain 
response. 
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304 Stainless Steel 


The problem of predicting low-cycle fatigue life at 
a geometrical notch has been studied by Krempl.? 
Circumferentially notched, cylindrical specimens and 
edge-notched flat plate specimens of 304 stainless steel 
were tested at room temperature and 550°F under 
fully reversed and zero-to-tension uniaxial loading 
conditions. Completely reversed bending tests of 
notched specimens were also performed at room 
temperature. In these tests, values of theoretical stress 
concentration factor, K;, ranged from 1.9 to 3.3. In all 
these tests, creep deformation was considered to be 
insignificant. Both the number of cycles to crack 
initiation (considered to be a crack 0.005 to 0.015 in. 
long) and to total failure of the specimen were 
determined. 

In the range of fatigue lives between 10? and 10° 
cycles, no shakedown to elastic action at the notch 
root was observed, and for the cylindrical specimens 
about 70% of the total life was involved in crack 
propagation.?* Knowing the applied nominal elastic 
stress amplitude and mean stress, the value of K;,, 
fatigue resistance of unnotched material, and the 
cyclic-stress—strain behavior of the material, a method 
was developed to predict fatigue life of the notched 
specimens that had been tested under load-controlled 
conditions. The method provided conservative predic- 
tions that may be of design use but did not provide 
accurate predictions. Additional experimental work at 
both room and elevated temperatures is required 
before a more accurate understanding of notched 
low-cycle fatigue behavior can be developed. 

Studies at GE-—NSP,@® at BCL,'%'7 at 
ANL,! 8?7 and at ANC?8"3! have all been conducted 
on material from one heat (No. 55697) of annealed 
(1 hr at 1950°F) 304 stainless steel; ANL has also 


Table 3 Comparison of Monotonic- and Cyclic-Stress—Strain Properties of Solution- 
Annealed Alloy 800 at a Strain Rate of 4 x 10°° sec! (From Jaske et al.'?) 





0.2% 











Max, , Ratio of 
total offset tensile ‘adiline 
i yield strength, ee Strength coefficient, Strain- 
strain ; monotonic 
ksi ; ksi hardening exponent 
range, Temp., yield 
% °F Monotonic Cyclic strength Monotonic Cyclic Monotonic Cyclic 
3.0 70 24.9 56.6 Leet 59.2 106 0.148 0.099 
3.0 800 IS DO.0 3.49 41,2 118 0.156 0.120 
3.0 1000 15.0 61.3 4.08 30.1 115 0.113 0.098 
3.0 1200 12.2 5355 4.38 27.9 92:1 0.123 0.087 
3.0 1400 1235 32,5 2.60 28.8 62.1 0.100 0.101 
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Fig. 12 Effect of irradiation on cyclic-stress—strain behavior 
of Alloy 800 at 1292°F and at a strain rate of 8 x 10% sec’. 
(From Brinkman et al.°!) 


conducted a few tests at 1100°F on specimens from a 
similar heat (No. 9T2796) of this alloy.?® In the work 
at GE—NSP, specimens were machined to final con- 
figuration and then reannealed (', hr at 1998°F).° 
Specimens at BCL were tested in the as-received 
annealed condition.'*"!7 Material tested at ANL was 
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in the same annealed condition used by GE—NSP as 
well as in a stress-relieved condition (‘, hr at 1400°F) 
and in an aged condition (159 and 1000 hr at test 
temperature prior to testing).25-?7 At ANC, material 
has been tested in an annealed (', hr at 1967°F) 
condition, in several aged conditions, and in irradiated 
conditions.?*3! 

Combining results from all four laboratories, 
fatigue curves can be defined for 304 stainless steel at 
70, 800, 1000, 1200, and 1500°F and at a strain rate 
of 4x 10° sec’. There is a general trend of decreased 
fatigue resistance with increased temperature. How- 
ever, for the annealed material, there is only a slight 
difference in fatigue strength between 1000 and 
1200°F (Ref. 14). Work at ANL?* showed that aged 
material (1000 hr at 1050°F) and_ stress-relieved 
material tested at 1050°F had significantly better 
fatigue strength than annealed material tested at 
1050°F. However, at 1200°F all three conditions had 
essentially the same fatigue resistance. This effect of 
heat treatment was related to precipitation of M,3C¢. 
carbides. Higher elevated-temperature fatigue resistance 
was related to more advanced carbide precipitation. 
Thus, for longer testing times and higher temperatures, 
the influence of initial heat treatment is expected to 
have a small effect on fatigue life because precipitation 
will have time to stabilize relatively early in the test. 

Equation | was again used to define best-fit fatigue 
curves as summarized in Table 4. Data on solution- 


Table 4 Summary of Best-Fit Regression Coefficients for 
304 Stainless-Steel Fatigue Curves 





Best-fit value for 


BEE: References from 
indicated constant 








Temp., Material Type of which data 
°F condition strain range Ay A, were taken 
70 Annealed dep —0.608 —0.403 16 
Aée —1.13 —0.303 
800 Stress relieved A€ép —0.422 —0.451 7 
Aée —1,18 —0.191 
1000 Annealed Ae€p —0.0429 —0.610 14 
Aée —1.02 —0.208 
1100 Annealed Aen —0.620 —0.467 25 
Aée —1.18 —0.176 
1100 Aged 1000 hr Aep —0.541 —0.441 25 
at 1100°F 
Meg 1,86 —0.177 
1200 Annealed and Aep +0.113 -0.681 7,14 
stress relieved 
Aée -1,88 —0.188 
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annealed material at 1100°F are for both heats that 
were tested by ANL, and data for aged material are for 
the second heat (No. 9T2796).?* Berling and Slot? 
have shown that decreasing the strain rate from 
4x 10° sec! to4 x10 and 4 x 10° sec! caused a 
corresponding reduction in cyclic- fatigue life at both 
1202 and 1501°F. (See Fig. 13.) 
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Fig. 13 Effect of strain rate on fatigue life of 304 stainless 
steel at 1202 and 1501°F, (From Berling and Slot.”) 


Irradiation of 304 stainless steel to fluences be- 
tween 1.5 and 3.0x 107? neutrons/em? (£>0.1 
MeV) at 1292 to 1382°F caused about a factor of 2 
reduction in fatigue life?! at 1292°F. (See Fig. 14.) 
The modified universal-slopes equation*®?’?? predicted 
a Similar reduction in fatigue life based upon changes in 
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Fig. 14 Effect of irradiation on the fatigue life of 304 
stainless steel at 1292°F and at a strain rate of 8 x 10% sec’. 
(From Brinkman et al.?') 


irradiated tensile properties and gave a conservative 
prediction of actual fatigue life. 


At both 1000 and 1200°F, tension hold times 
caused a significant reduction in fatigue life compared 
to results for continuous cycling. At 1200°F, compres- 
sion and combined tension—compression hold times 
did not significantly affect cyclic-fatigue life. Results at 
1200°F were for annealed specimens tested at GE— 
NSP,® and those at 1000°F were for annealed speci- 
mens tested at BCL.°® As was true in the case of Alloy 
800, the tension hold-time data can be represented by 
a series of curves for each constant length of hold time. 
(Refer to Fig. 15 and Table 5). For hold times up to 
5 hr at 1000°F, cyclic-fatigue resistance was reduced as 
length of hold time increased [Fig. 15(@)]; for hold 
times up to }, hr, the same was true at 1200°F 
[Fig. 15(b)]. However, at 1200°F there was a satura- 
tion in the hold-time effect for hold times from ', to 
10 hr, where there was little difference in fatigue life as 
length of hold time increased. This saturation is 
probably related to the enhanced carbide precipitation 
at the grain boundaries discussed earlier. 


For this alloy, Eq. 2 was also found to describe 
cyclically stable relaxation behavior both at 1000°F 
(Ref. 36) and at 1200°F (Ref. 37). Equation 3 was 
used to calculate creep—fatigue damage at 1000°F 
(Ref, 36) and at 1200°F (Ref. 38), and the results are 
summarized in Fig. 16. The damage values fall near or 
well below the D=1 curve, and the trend of values is 
different for different strain ranges. (No such trend was 
noted for Alloy 800.) Fatigue damage tended to 
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Fig. 15 Effect of hold time on fatigue life of 304 stainless steel at a strain rate of 4 x 10° sec’, 
(From Conway? and Jaske et al,*°) 


Table 5 Slope and Intercept Values That Define 
Tension Hold-Time Fatigue Curves for 
Annealed 304 Stainless Steel 








Length of 
tension Slope of Ac; References from 
Temp., Intercept hold time, vs. Ne which data 
°F at Nr=1 min curve were taken 
1000 0.83 0 —0.56 14 
10 —0.61 36 
30 —0.67 36 
60 —0.69 36 
300 —0.81 36 
1200 0.62 0 —0.54 7, 14, 37 
1 —0.59 8, 37 
10 —0.65 8, 37 
30 to 60 —0.69 8,37 
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Fig. 16 Creep—fatigue damage interaction for tension hold- 
time tests of annealed 304 stainless steel. (From Jaske et al.°° 
and Campbell.**) 


decrease toward a limit as the length of hold time 
increased, whereas the creep damage tended to first 
increase slightly and then decrease rapidly. Total- 
damage values decreased as a function of hold time as 
illustrated in Fig.17. Type 304 stainless steel was 
irradiated to fluences of 1.65 to 5.02 x 10?! neutrons/ 
cm? (E >0.1 MeV) at 842°F and tested with tensile 
hold times at Ae, =0.01 and at 1100°F (Ref. 30). 
Compared with control specimens that were aged at 
752°F for 8760 hr and at 770°F for 3096 hr, irradia- 
tion reduced the tension hold-time fatigue life by a 
factor of approximately 4. Total damage, as calculated 
by Eq.3, was between 1 and 3 for the control 
specimens with tensile hold times from 0.01 to 1.0 hr. 
However, the total damage was between 0.17 and 0.72 
for irradiated specimens with hold times from 0.1 to 
5.0 hr. 

Incremental step tests were used to develop the 
monotonic- and cyclic-stress—strain curves shown in 
Fig. 18 (Refs. 12 and 15). At high strain amplitudes, 
results of the incremental step tests correlate fairly well 
with the constant-amplitude data points. However, at 
low strain amplitudes, the step-test curves reflect the 
fact that the material was previously hardened at a 
higher level and is thus harder than observed in 
constant-amplitude data. This same “memory” be- 
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havior of the material was shown by less cyclic 
hardening for specimens from step tests with lower 
maximum strain amplitudes. (Compare curves for 
Ae;/2 = 0.015 max. to those for Ae;/2 = 0.005 max. in 
Figs. 18(c) and (d).) It was also shown that decreasing 
the strain rate from 4 x 10°° to 4 x 107° sec! caused 
more cyclic hardening at 1000°F and less at 1200°F. 
These results imply that the hardening mechanism is 
time dependent as well as cycle dependent. The 
time-dependent hardening is probably related to the 
formation of carbides discussed earlier. Values of the 
ratio of cyclic to monotonic yield strength (Table 6) 
show that this alloy hardened by factors from 1.9 to 
3.0, with maximum cyclic hardening at 1000°F. 
Equation 4 was again found to well represent both 
monotonic- and cyclic-stress—strain curves, and the 
appropriate constants are listed in Table 6. Irradiation 
of 304 stainless steel to 3.5 x 107! neutrons/cm? 
(E >0.1 MeV) at 842°F was shown to only slightly 
increase the amount of cyclic hardening at 1100°F 
(Ref. 30). 


316 Stainless Steel 


Material from one heat (No. 65608) of annealed 
316 stainless steel has been used in low-cycle fatigue 
programs at GE-—NSP,® BCL,'%!7 ANL,?° and 
ANC.?®3! Specimens used at GE—NSP were rean- 
nealed (/, hr at 1958°F) and then stress relieved (1 hr 





102 Parry =z. a a 
©, 1200°F, Ae,=0020 } ‘ 

A. 1200°F, Ae=0005 f CAmPbE 

@, 1000°F, Ac,=0.022 

B, 1000°F, Ac,=O.010 a 

A, 1000°F, Ae,=0004 











8 iL 4 
rites : 
4 i Trend at Ae=O0.020 at 1200°F 

= | Trend at Ae=O005 at 1200°F 

if a Trend at Ae,=0.022 at {(000°F 

<a D=1 

E i? x 

lo 


2 















ee? oe 
10° 10! 107 10° 
LENGTH OF TENSION HOLD TIME (t,), min 


Fig. 17 Effect of length of tension hold time on total damage 
of annealed 304 stainless steel. (From Jaske et al.°° and 
Campbell.?* ) 
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Fig. 18 Comparison of monotonic- and cyclic-stress—strain curves for annealed 304 stainless steel, 


(From Jaske et al,'?°!5) 


at 1400°F) after final machining.® This same treatment 
was used for some of the specimens tested at 1050°F 
by ANL.?3°?5 ANL has also used specimens in the 
annealed condition at both 1050 and 1202°F and 
specimens that were aged (1000 hr at 1050°F) in tests 
at 1050°F (Ref. 25). All specimens tested at BCL were 
either in the annealed condition or were aged for 
1000 hr at either 1050 or 1200°F prior to test- 
ing.'5“!7 Annealed, aged, and irradiated materials have 
been used in the experimental work at ANC.?*3! 

For constant-amplitude, continuous cycling of 316 
stainless steel, fatigue curves are available at 70, 800, 


900, 1050, 1200, and 1500°F. With increasing tem- 
perature, there is a general trend of decreasing low- 
cycle fatigue resistance. At both 1050 and 1200°F, 
there was no significant difference in fatigue strength 
between stress-relieved specimens tested at ANL?? and 
annealed specimens tested at BCL.'® Material that was 
aged at test temperature for 1000 hr before testing at 
1050 and 1200°F had slightly better fatigue strength 
than the other two heat-treatment conditions at strain 
ranges above 0.005 (Refs. 16 and 17). At strain ranges 
below 0.005, the aging had no significant effect on 
fatigue life. Results of tests at 1100°F showed a similar 
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Table 6 Comparison of Monotonic- and Cyclic-Stress—Strain Properties for Annealed 
304 Stainless Steel (From Jaske et al.') 














0.2% , 

—. offset tensile _ = Strength Strain- 

total é cyclic to ips ; 

title Strain yield strength, deimaiiaaie coefficient, hardening 

range, Temp., rate, ksi yield kel ee 
%o BR i oO Monotonic Cyclic strength Monotonic Cyclic Monotonic Cyclic 
3.0 70 4x103 31.4 66.3 21 RS a 0.152 0.163 
3.0 800 4x10? 17.5 48.3 2.76 Sted 98.9 0.126 0.116 
3.0 1000 4x10? 14.6 43.5 2.98 43.6 80.6 0.164 0.097 
1.0 1000 4x10? 14.6 34.3 2:39 86.7 82.6 0.293 0.143 
1.0 1000 4x10% 14.6 40.0 2.74 87.0 98.2 0.288 0.141 
3.0 1200 4x10? 14.2 36.6 2.58 41.8 65.2 0.172 0.101 
1.0 1200 4x10? 14.2 32.4 2.28 81.8 68.4 0.182 0.120 
1.0 1200 4x10% 1352 24.6 1.86 40.4 44.8 0.183 0.155 





influence of aging (either 3096 hr at 1112°F or 8760 hr 
at 1382°F on low-cycle fatigue resistance of 316 
stainless steel.2° The effect of strain rate has been 
reported by Berling and Slot” and is similar to that 
discussed earlier for 304 stainless steel. Best-fit fatigue 
curves were defined using Eq. 1, and the results are 
listed in Table 7. 

When 316 stainless steel was irradiated to fluences 
between 1.3 and 2.4 x 107? neutrons/em? (E>0.1 
MeV) at 1292 to 1382°F, fatigue life at 1292°F was 
reduced by a factor of approximately 2, which was 
similar to comparable data for 304 stainless steel.*' 


(See Fig. 19.) Just as for the 304 alloy, the modified 
universal-slopes equation??’? predicted a similar 
amount of reduction in fatigue life because of the 
difference between tensile properties of the irradiated 
and unirradiated material, and the prediction of actual 
fatigue life was conservative. 

As was true for both Alloy 800 and 304 stainless 
steel, experimental work has shown that tension hold 
times detrimentally affect the cyclic-fatigue life of 316 
stainless steel®*!®*!7°3° at temperatures from 1050 to 
1200°F. In the same temperature range, limited tests at 
Ae; = 0.02 have shown that compression and combined 


Table 7 Summary of Best-Fit Regression Coefficients 
for 316 Stainless-Steel Fatigue Curves 





Best-Fit value of 








Type of fies References from 
Temp., Material strain a which data 
°F condition range A, A, were taken 
70 Annealed Aep —0.176 —0.499 17 
Aée —1.75 —0.179 
800 Annealed Aep —0.374  —-0.484 tyhd 
Aée —1.62 —0.219 
900 Annealed A€p —0.396 -0.495 17 
Aée —1.74 —0.196 
1050 Annealed Aep —0.552 —0.497 17 
Aée —1.76 —0.186 
1050 Aged for 1000 hr Aep 0.209 —-0.727 17 
at 1050°F hee —1.97 —0.120 
1200 Annealed Aep —0.264  -0.560 4547 
Aée -1.91 —0.153 
1200 Aged for 1000 hr Aep 0.0871  -—0.680 17 
at 1200°F Aée -1.88 —0.180 
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Fig. 19 Effect of irradiation on the fatigue life at 316 stainless 
steel at 1292°F and at a strain rate of 8 x 10% sec'. (From 
Brinkman et al.?!) 


compression and tension hold times have no significant 
effect on cyclic-fatigue life.'7 The effect of tension 
hold time on fatigue life is illustrated in Fig. 20. For 
annealed material at 1050 and 1100°F [ Figs. 20(a) and 
(b)], fatigue resistance decreased as the length of ten- 
sile hold time increased. (See Table 8.) However, at 
1200°F, a saturation in this kind of behavior trend was 
noted, and similar fatigue lives were obtained for ten- 
sion hold times from 0.1 to 1.0 hr. This saturation be- 
havior is similar to that described previously for 304 
stainless steel at 1200°F. 

Aging of specimens for at least 1000 hr at test 
temperature prior to testing provided additional insight 
on the hold-time fatigue behavior of this alloy. At 
1050°F [Fig. 20(@)], aging had little influence on the 
cyclic-fatigue life of tension hold-time specimens. At 
1100°F [Fig. 20(b)], aging improved the hold-time 
fatigue life to a noticeable degree, and, at 1200°F 
[Fig. 20(c)], aging produced a very large increase in 
hold-time fatigue resistance. The extremely beneficial 
influence of aging at the highest of these three tempera- 
tures is probably related to the type of carbide pre- 
cipitation discussed previously for the 304 alloy. 
Further study is required to see if this trend is limited 
to the 1200°F temperature region or if it is carried on 
to higher temperatures. 

Both ANC?° and BCL!” have used Eq. 2 to 
describe the cyclically stable stress-relaxation response 
to this alloy and have computed creep—fatigue stress- 
relaxation response of this alloy and have computed 
creep—fatigue damage interaction by means of Eq. 3. 
As shown in Fig. 21, damage values for the annealed 
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material fell somewhat below the D=1 curve. The 
dispersion of points was not a definite function either 
of strain range, as was true in the case of 304 stainless 
steel, or of hold time (see Fig. 22), as in the case of 
both Alloy 800 and 304 stainless steel, but did reflect 
normal scatter in such experimental results. At 
1050°F, damage values for aged material were about 
the same as those for annealed material. However, at 
1100 and 1200°F, damage values for aged material fell 
close to or significantly above the D=1 curve. The 
hold-time fatigue lives of specimens that were irra- 
diated to fluences of 1.7 to 28.6 x 10?° neutrons/cm? 
(E >0.1 MeV) at 842°F were about a factor of 2 to 4 
below those of control specimens that were aged for 
8760 hr at 752°F. In terms of total damage calculated 
by means of Eq. 3, similar results were obtained for 
both the irradiated and control specimens because the 
damage calculation did account for the effect of 
irradiation on fatigue life. For hold times from 0.01 to 
1.0 hr, damage values between 0.51 and 0.85 were 
computed for the control specimens, and damage 
values between 0.28 and 0.85 were computed for the 
irradiated specimens. 

Cyclic-stress—strain curves from incremental-step 
tests of 316 stainless steel’? are presented in Fig. 23. 
The overall data trends for this alloy are similar to 
those discussed for 304 stainless steel. High-level strain 
cycles caused the stable stress on subsequent low-level 
cycles to be higher than for constant-amplitude load- 
ing. This phenomenon was more evident at room 
temperature [Fig.23(a)] than it was at elevated 
temperatures. The cyclic-stress—strain “memory” of 
the material resulted in less cyclic hardening for tests 
that were conducted at lower maximum-strain ampli- 
tudes. 

The complexity of cyclic-hardening behavior is 
further pointed out by the influence of aging and strain 
rate on cyclic hardening at 1050 and 1200°F. At 
1050°F [Fig. 23(c)], the aged material showed slightly 
less cyclic hardening than the annealed material, 
whereas no significant difference in monotonic-stress— 
strain response was noted. Decreasing the strain rate 
from 4 x 10° to 4.x 10° sec? resulted in both more 
monotonic and more cyclic hardening at this tempera- 
ture. At 1200°F [Fig. 23(d)], aging significantly re- 
duced the amount of cyclic hardening and had little 
effect on the monotonic response. When strain rate was 
decreased from 4x 10° to 4x 10° sec’, more 
monotonic hardening was exhibited, but significantly 
less cyclic hardening was evident. 

As shown by the ratios of cyclic to monotonic 
yield strength in Table 9, this alloy cyclically hardened 
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Fig, 20 Effect of hold time on fatigue life of 316 stainless steel at a strain rate of 4 x 10°? sec’, 
(Parts a and c from Jaske et al.!7 and Part b from Brinkman et al.°°) 


by factors between 1.3 and 2.6 with a maximum 
amount of hardening at 1050°F. Both monotonic- and 
cyclic-stress—strain curves were well represented by 
Eq.4 using the appropriate constants in Table 9. 
Irradiation of 316 stainless steel to 1.2 x 10?! 
neutrons/cm? (E >0.1 MeV) at 842°F caused only a 
slight increase in cyclic hardening at 1100°F (Ref. 30). 
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SUMMARY 


Data from axially loaded, strain-controlled tests of 
unnotched specimens are being used to evaluate the 
low-cycle fatigue behavior of reactor structural 
materials and to provide creep—fatigue design informa- 
tion on these materials. For continuous cycling (with- 
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Table 8 Slope and Intercept Values That Define 
Tension Hold-Time Fatigue Curves for 
Annealed 316 Stainless Steel 





Length of Slopeof Ac, References from 








Temp., Intercept tension hold vs. which data 
°F at Nr=1 time, hr Nf curve were taken 
1050 0.214 0 —0.40 17 
0.1 —0.50 
1.0 —0.55 
5.0 —0.63 
1100 0.284 0 —0.44 30 
0.01 —0.48 
0.10 —0.56 
0.50 —0.63 
1.0 —0.72 
Sun —0.83 
10.0 —0.97 
1200 0.293 0 —0.45 17 
0.1 to 1.0 —0.60 
out hold times), fatigue curves have been developed for steel have similar low-cycle fatigue strength, whereas 
Alloy 800 and 304 and 316 stainless steels over 316 stainless steel has slightly lower fatigue strength 
temperatures from 70 to 1500°F at a strain rate of than the 304 alloy.?® Irradiation reduced the fatigue 
4x 10° sec'. Decreasing the strain rate down to life of the stainless steels by a factor of about 2 and of 
4x 10° sec’ caused significant reduction in the the Alloy 800 by factors up to 35 at 1292°F (Ref. 31). 
low-cycle fatigue resistance of these alloys at tempera- Tension hold times had a detrimental effect on the 
tures of 1200 and 1500°F (Ref. 7). In the 1100 to low-cycle fatigue resistance of all three alloys at 
1200°F temperature range, Alloy 800 and 304 stainless temperatures of 1000°F and greater. In contrast, 


compression or combined compression—tension hold 
times had little or no effect on low-cycle fatigue 
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Fig. 21 Creep-fatigue damage interaction for tension hold- Fig. 22 Effect of length of tension hold time on total damage 
time tests of 316 stainless steel. (Results at 1050 and 1200°F of 316 stainless steel. (Results at 1050 and 1200°F from Jaske 
from Jaske et al.'7 and at 1100°F from Brinkman et al.°°) et al.!7 and at 1100°F from Brinkman et al.*°) 
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Fig. 23 Comparison of monotonic- and cyclic-stress—strain curves for 316 stainless steel. (From Jaske 
et al,'7) 
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Table 9 Comparison of Monotonic- and Cyclic-Stress—Strain Properties 
of 316 Stainless Steel (from Jaske et al.' 7) 





Max, Ratio of 











0.2% offset ; Strength Strain- 
total : cyclic to fh : 
; : yield strength, ; coefficient, hardening 
Strain strain : monotonic : 
ksi ‘ ksi exponent 

Temp., rate, range, yield 
°F sec % Monotonic Cyclic strength Monotonic Cyclic Monotonic Cyclic 
70 4x10" 5.78 39.4 80.6 2.05 86.1 146 0.095 0.099 
70 4x10° 2.90 39.4 69.9 Lk 84.2 133 0.086 0.109 
800 4x10% 2.93 26.8 63.1 235 47.4 118 0.102 0.104 
800 4x10° 1.00 26.8 39.8 1.48 62.3 104 0.140 0.154 
1050 4x10 2.92 22.3 58.9 2.64 40.9 102 0.100 0.089 
1050 4x 10° 1.00 223 45.2 2.03 55.8 91.1 0.142 0.113 
1050 4x10° 0.99 24.1 $3.1 23 46.0 100 0.106 0.101 
1050* 4x10% 2.87 222) 56.8 255 46.2 92.3 0.122 0.078 
1050* 4x10°% 0.98 22.3 43.7 1.96 46.3 Bg 0.119 0.151 
1050* 4x10% 1.00 24.1 49.2 2.04 43.2 96.4 0.096 0.109 
1200 4x10° 2.90 20.3 50.3 2.48 48.3 99.7 0.142 0.115 
1200 4x10° 0.97 20.3 40.5 1.99 47.8 71.1 0.143 0.090 
1200 4x10° 0.98 22.9 Be 1.54 42.8 67.3 0.102 0.106 
1200t 4x10°% 2.90 20.3 41.5 2.04 54.6 Ta 0.159 0.094 
1200; 4x10% 0.98 20.3 34.2 1.68 49.1 70.9 0.112 0.116 
1200t  4x105 1.00 21.6 28.8 1.33 52.9 48.0 0.147 0.081 





*Aged 1000 hr at 1050°F. 
+Aged 1000 hr at 1200°F. 


resistance. In terms of hold-time fatigue resistance of 
the annealed material, 304 stainless steel was superior 
to 316 stainless steel, which was in turn superior to the 
Alloy 800. Since aging had a large beneficial influence 
on hold-time fatigue behavior of both 316 stainless 
steel and Alloy 800 at 1200°F (Ref. 17), a direct 
comparison of the three alloys cannot be adequately 
made until more data on the influence of aging are 
available. The effect of aging is related to carbide 
precipitation,?* and the beneficial influence of aging 
corresponds directly with the amount of cyclic harden- 
ing observed (i.e., less hardening was observed when 
aging was beneficial).!7 Because experimental results 
indicate that aging is probably beneficial when it has an 
effect on elevated-temperature creep—fatigue re- 
sistance, information developed on annealed material 
will provide conservative estimates for design use. 
However, it should be pointed out that, although aging 
may inhibit high-temperature damage mechanisms, it 
may enhance fatigue damage at lower temperatures. It 
has been shown that heat treatment of 316 stainless 
steel to produce 2 to 3 vol.% (CrFe)3C, carbides will 
reduce the room-temperature low-cycle fatigue re- 
sistance as compared with  solution-annealed 
material.?® Thus the influence of combined thermal— 


mechanical (rather than just isothermal) cycling on 
fatigue resistance of these alloys should be investigated 
in future work, 

Use of the linear-life-fraction damage rule to 
evaluate creep—fatigue damage interaction gave dif- 
ferent results for each of the three alloys. Calculated 
damages values varied as a function of both strain range 
and hold time for 304 stainless steel, varied as a 
function of hold time but not of strain range for Alloy 
800, and did not vary as a function of either strain 
range or hold time for 316 stainless steel. Such 
different trends indicate that relating creep damage 
done during cyclic relaxation to monotonic-creep- 
rupture behavior does not explain the time-dependent 
damage done in hold-time fatigue tests. 

The incremental step test provided a useful method 
of investigating the cyclic-stress—strain behavior of 
these alloys. The cyclic-stress—strain response of all 
three materials at low strain levels was definitely 
influenced by prior cyclic hardening at high strain 
levels. Such behavior was more pronounced for the 
stainless steels than it was for the Alloy 800. This 
influence of variable-amplitude loading is an important 
consideration in design applications where materials are 
subjected to spectra of loads rather than constant- 
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amplitude loading. All three alloys cyclically hardened 
by a significant amount from the annealed condition. 
Alloy 800 hardened by factors from 2.3 to 4.4 witha 
maximum amount of 1200°F, 304 stainless steel 
hardened by factors from 1.9 to 3.0 with a maximum 
amount at 1000°F, and 316 stainless steel hardened by 
factors from 1.3 to 2.6 with a maximum amount at 
1050°F. The cyclic-stress—strain response of both 
stainless steels was influenced by strain rate. (The 
effect of strain rate on Alloy 800 was not examined.) 
At 1000 to 1050°F, decreasing the strain rate caused 
increased cyclic hardening, and, at 1200°F, it caused 
decreased cyclic hardening. 
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ANS GUIDELINES FOR PREPARING CRITICAL REVIEWS 
FOR REACTOR TECHNOLOGY 


One interpretation of the words “Critical Review’ emphasizes the word critical. That type of 
article would have as its purpose the discussion of a single subject, which has been uncertain and 
perhaps controversial, in considerable depth. An example might be a critical review of the values of 
alpha for plutonium, assessing all the work done and arriving at a best current estimate. This paper 
is primarily addressed to specialists on the particular subject, and serves as an authoritative source 
for the information they use. 

In the second interpretation of a Critical Review, the emphasis is on the word review, in the 
sense of survey. Such a paper would be broader and probably more descriptive in scope. A paper 
on ‘Solubility of Metallic Elements in Liquid Sodium” is an example of this category. Such a 
paper would be addressed to a much broader group of readers and written in a fashion that would 
be of interest to the majority of reactor technologists, |i would provide enough of an introduction 
that the specialist from another field could immediately appreciate why the subject is important. 

The criteria for these two types of articles may not be very different: 

1. The paper should be based on a thorough coverage of relevant work on tie subject from 
many sources, A review based on the work of only one individual or group is better suited for 
publication in one of the regular society journals. The review should be critically selective, 
reporting only the most valid results and indicating why some prior results have a questionable 
status. The review should call attention to significant gaps where more work is required in the 
subject of the article. 

2. The paper must include an evaluation, a critical appraisal, of the field being covered. It 
should be more than merely a survey of the work being done in the field. 

3. The paper should be timely, on a subject of active current interest. 

4. The scope of articles acceptable as Critical Reviews includes all the subject areas identified 
for Reactor Technology (Economics, Physics, Mechanics, Construction, Fuel Elements, Fuel 
Cycles, Fluid and Thermal Technology, Fuel Processing, Components, Operating Performance), as 
well as Materials (including Source and Special Nuclear), Environmental Effects, and Effluent 
Management, Not desired are (1) reports of original research proposed for first publication and 
(2) review articles directed toward the specialist in the field of nuclear safety (which is covered by 
the AEC’s bimonthly review Nuclear Safety). 

5. The paper should be organized so that it is of immediate practical use to the readers. Such 
Organization requires: attention to consistent use of units, presentation of important data in 
summary curves or tables, and a fully adequate bibliography to the original literature. 


All papers are reviewed, at the proposal and the draft stage, by several referees expert in the 
subject field. Referees may require revisions, expansions, or deletions to make the paper acceptable 
for publication. Details on size, honorarium, style, etc., can be obtained from ANS. 
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Flux Synthesis Methods 
in Reactor Physics 


By Weston M. Stacey, Jr.* 


Abstract: Use of flux synthesis methods to compute the 
neutron distribution and other important parameters in nuclear 
reactors is reviewed, Spatial, spectral, and angular synthesis 
methods are discussed, Examples are described which illustrate 
common methods of selecting trial functions and which 
represent the current state of the art. Recent advances in the 
underlying theory and in the development of specific applica- 
tions are discussed, experience with the different synthesis 
methods is assessed, and promising areas for further develop- 
ment activity are suggested. 


Synthesis, which derives from the Greek, means 
literally the action of putting together. Webster defines 
it as the “composition or combination of parts or 
elements so as to form a whole.” Thus, in a literal 
sense, most of the methods that comprise the 
repertoire of the reactor physicist may be classified as 
synthesis methods. Such a classification, although 
serving to unify a minuscule portion of the universe, 
has little pragmatic advantage. To most reactor 
physicists the term “flux synthesis” conjures up a 
vaguely delineated class of methods having to do with 
the combining of simpler or lower dimensional solu- 
tions to obtain a more complex or higher dimensional 
solution, 

Meyer! seems to be the first to have published a 
synthesis technique for constructing multidimensional 
flux solutions from lower dimensional calculations. His 
idea was to calculate a different two-dimensional 
planar flux solution for each of several axial zones. 
These planar solutions were then used to obtain 
planar-averaged cross sections and transverse bucklings 
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for the corresponding axial zones, and one or more 
axial shapes were calculated corresponding to the 
planar regions over which the cross sections and 
transverse bucklings were computed. Finally the planar 
and axial shapes were multiplied together to yield a 
three-dimensional flux distribution. Wachspress et al.” 
improved on this by using the planar solutions to 
define coupling coefficients between contiguous planar 
regions (channels), rather than transverse bucklings. 
Thus, in their multichannel synthesis technique, the 
axial-shape calculations (for the axial shapes corre- 
sponding to the different planar regions or channels) 
are coupled rather than independent as in Meyer’s 
method. 

Selengut® precipitated the subsequent development 
of synthesis methods by providing a_ theoretical 
foundation based on variational theory and by intro- 
ducing the idea of “blending” or combining two or 
more lower dimensional solutions. Calame and 
Federighi* were the first to exploit Selengut’s ideas. 
They derived variationally a method to synthesize 
“hard” and “soft” Maxwellian spectra to obtain 
spatially dependent thermal-neutron spectra. Their 
method was originally referred to as the overlapping 
group method but has more recently come to be 
designated the spectral synthesis, or space—energy 
synthesis, method. 

Kaplan* employed Selengut’s idea of blending to 
extend Meyer’s technique. In Kaplan’s method the 
planar solutions calculated for the different axial zones 
were all used over the entire axial domain, rather than 
each planar solution being used only in the correspond- 
ing axial zone as in Meyer’s method. Axial-shape 
functions corresponding to the different planar solu- 
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tions were determined by solving a coupled set of 
axial-shape equations. In Kaplan’s method a single axial- 
shape function was associated with each planar solu- 
tion, and all the planar solutions were used in every 
axial zone, which has given rise to the designation 
single-channel continuous spatial synthesis. Kaplan also 
noted that the synthesis equations could be derived on 
the basis of the weighted residual method, as well as 
variationally, and that the equations obtained from the 
two derivations were essentially equivalent. 

Pomraning and Clark® appear to be the first to 
have exploited Selengut’s ideas with respect to 
synthesizing the angular distribution of the neutron 
flux. They demonstrated that the synthesis equations 
were, in fact, the Py equations when the Legendre 
polynomials were used to expand the angular flux and 
that the use of a different set of polynomials led, via 
the synthesis formalism, to an improved form of 
diffusion theory. 

An idea of Selengut’s’ was also the basis for the 
next advance in the development of synthesis methods. 
He proposed a variational principle that would admit 
discontinuous trial solutions. Wachspress and Becker® 
exploited this idea to combine the blending formalism 
of Kaplan and the earlier multichannel idea of 
Wachspress et al. to obtain a multichannel synthesis 
formalism in which the planar solutions were blended 
with axial shapes calculated for each planar region or 
channel. These axial shapes satisfied coupled equations. 
Yasinsky and Kaplan? also employed Selengut’s varia- 
tional principle to derive an axially discontinuous 
extension of Kaplan’s single-channel synthesis in which 
different sets of planar shapes were blended in the 
different axial zones. 

Subsequent to these pioneering papers, the litera- 
ture of the field has been swollen with reports of 
development and application of synthesis methods and 
extensions of the supporting theory. It is the purpose 
of this article to review this work critically and 
selectively, with emphasis on that work which provides 
a measure of the current state of the art. The review 
also touches on several areas in which essentially 
synthesis-like procedures were employed in the deriva- 
tion of approximations. 

There is no generally accepted classification of 
synthesis methods or nomenclature for the different 
approximations. However, to facilitate the organization 
of the article, it is convenient to establish some 
conventions before proceeding. To this end, those 
methods which most logically may be characterized as 
the combining of spatial shapes, energy spectra, or 
angular functions will be referred to as spatial, spectral, 
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and angular synthesis, respectively. Some typical 
synthesis approximations are so classified below. 


Spatial Synthesis: 
N 
$(x,y,2,E,Q,1) = L VnlXy)Pn(Z,£,2,0) (1) 


or 
N 
L Vn(X,¥,Z) Pn (E,Q,0) (2) 
or 
N 
Xx Vn(x,y.E)pn(z,2,0) (3) 


Spectral Synthesis: 
N 
$(x,y,2,E,Q,1) = Lo Yn(E)on(x,¥,2,2,0) (4) 
n= 


Angular Synthesis: 
N 
(x,y,2,E,Q,t) = Le Vn(2Q)Pn(x,y,2,E,t) (5) 
R= 


where W,, = known expansion function 
Pn = unknown expansion coefficient 
x,y,z = spatial coordinates 
FE = energy 
Q2 = direction 
t = time 


It should be emphasized that the synthesis approxi- 
mations are generally based on discrete approximations 
(e.g., multigroup treatment of energy, finite-difference 
treatment of space, and Py treatment of angle), rather 
than on the integro-differential Boltzmann equation; 
however, it is notationally convenient to indicate 
continuous variables above. 

In the derivation of synthesis approximations, it is 
necessary to select weighting functions as well as trial 
solutions (or trial functions as they are frequently 
called). The synthesis equations, in essence, are require- 
ments that the combining coefficients, with which the 
trial functions are combined, are selected such that the 
synthesized solution satisfies the neutron balance 
equation in a weighted integral sense, for as many 
different weighting functions as there are trial func- 
tions. Three different types of weighting functions are 
common: (1) Galerkin weighting, in which the weight- 
ing functions are identical with the trial functions; 
(2) adjoint weighting, in which the weighting functions 
are trial adjoint solutions; and (3) subdomain weight- 
ing, in which the weighting functions are unity over a 
subdomain and zero elsewhere. 
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SPATIAL SYNTHESIS 

The effort devoted to the development of spatial 
synthesis: far exceeds that expended-on the other 
categories of synthesis; thus it is not surprising that the 
greatest success has been realized in this area. Although 
the blending method of Kaplan has enjoyed the most 
extensive development, the nonblending method of 
Meyer has been successfully extended as well. The 
multichannel blending method, although promising, 
has never quite caught on, due to certain difficulties 
that may now be eliminated by a new formulation. The 
use of piecewise polynomials was introduced several 
years ago in the spatial synthesis context and has 
recently engendered considerable enthusiasm now that 
a more systematic development is taking place under 
the name of finite-element methods. Coarse-mesh 
rebalancing methods have recently, and correctly, 
come to be viewed as a variant of the spatial synthesis 
method. 


Nonblending Methods of Spatial Synthesis 


Nonblending methods, because of their relative 
simplicity, probably have been used more widely, if 
less extensively, than the more powerful blending 
methods. In the simplest of these methods, a single, 
but different, planar solution is used within each axial 
zone to obtain planar-averaged group constants for the 
calculation of an axial shape, which is then multiplied 
by the planar solutions to construct a_ three- 
dimensional solution, This method, which has probably 
been used at some time at every installation concerned 
with reactor analysis, can provide a useful approxima- 
tion in those situations in which the planar and axial 
flux distributions, to a good approximation, are 
separable within the different axial regions. Such 
situations become vanishingly few as calculational 
models become more realistic and reactor designs 
progress beyond their initial stages. 

Iteration between the axial- and _ planar-shape 
calculations is one means of improving the accuracy of 
the nonblending method. Travelli and Helm'° em- 
ployed an iterative nonblending synthesis method to 
calculate the spatial distribution of the sodium-void 
reactivity effect in the fast reactor critical facilities 
ZPR-6, Assemblies 2 and 3, which were large, single- 
zoned, reflected reactors. Both the real and adjoint 
fluxes were assumed separable in r and z coordinates 
over the whole reactor. One-dimensional 26-group 
calculations in r and z geometries were run through a 
rapidly converging iteration scheme by modifying the 
(uniform) transverse bucklings to obtain a pair of 
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calculations yielding identical kegg and with transverse 
bucklings which summed to the material buckling. 
Direct two-dimensional 26-group calculations were 
performed for comparison. The results are summarized 
in Table 1, and the spatially detailed comparison for 


Table 1 Synthesis!°® of Sodium-Void Reactivity 











Worth for ZPR-6 
Assembly 2 Assembly 3 
Synthe- Synthe- 
sis Direct sis Direct 
Total sodium void 
(10°? &k/k) —2.00 —2.03 —1.68 -1.71 
Central sodium void 
(10°? &k/k) —5.57 —5.68 -3.37 -3.46 
Ket 0.999 0.998 0.998 0.997 





Assembly 2 is shown in Fig. 1. The time required for 
this relatively simple synthesis method was only 1% of 
that for the direct two-dimensional solution, yet it 
yielded quite accurate results for this geometrically 
simple model. 

Hutchins et a employed an iterative extension 
of Meyer’s synthesis method to calculate the important 
properties of a 1000-MW(e) PuO,—UO, liquid-metal- 
cooled fast breeder reactor (LMFBR), using 16-group 
diffusion theory. Initial radial solutions were obtained 
through the core, followed by the calculation of 
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Fig. 1 Spatial distribution of sodium-void effect’ ° (—1 x 10°? 
Ak/k) in ZPR-6, Assembly 2, —, synthesis; - - -, direct 2D. 
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flux-adjoint weighted radial bucklings. Axial calcula- 
tions were performed in from one to four core zones 
using the bucklings from the initial radial solution. (Up 
to this point the method is identical to Meyer’s.) Next 
the axial regions were homogenized by flux-adjoint 
weighting with the axial solutions. A final radial 
solution was made, using the axially homogenized cross 
sections, and combined with the axial shapes to 
construct the two-dimensional flux. This synthesis 
solution was a factor of 20 faster than a direct 
two-dimensional solution of the same problem. 

Results of the synthesis and direct solutions are 
compared in Table 2. The distribution of power by 


Table 2 Results for a 1000-MW(e) Fast Reactor Using 
Synthesis and Conventional Two-Dimensional 
Solutions! ! 








Synthesis Conventional 
solution solution 
kee 0.9999 1.0002 
Total breeding ratio 1.338 1333 
Internal conversion ratios 
Core zone 1 0.9170 0.9172 
Core zone 2 0.7109 0.7112 
Percent power by region 
Core zone 1 50.4 50.7 
Core zone 2 39.2 39.0 
Radial blanket 2.81 2S 
Axial blanket TAS 7.06 
Corner blanket 0.46 0.51 
Maximum sodium loss 
reactivity, Ak/k +0.0235 +0.0221 





regions was accurate to 1% in the core and 10% in the 
corner blanket. Detailed power distributions in the 
core were accurate to with 5%. The maximum sodium 
loss coefficient was accurate to within 10%, and the 
Doppler coefficient was computed more accurately. 


Blending Methods of Spatial Synthesis 


Blending methods were formulated by Kaplan.* 
Subsequent theoretical developments have extended 
the methodology to time-dependent problems'? to 
accommodate the use of different sets of planar 
solutions in different axial regions’?’'*’7’ and to 
accommodate the use of different sets of planar 
solutions during different intervals of time.'®*'? The 
accuracy of the blending methods has been tested by a 
relatively large number of numerical studies performed 
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19°12 .20-24 25-27 


for both therma and fast reactor 
models. These tests have been generally successful, 
except for certain anomalies that will be discussed in a 
subsequent section, and have demonstrated the power 
of the synthesis blending methods. 

An example presented by Yasinsky and Kaplan? 
serves to represent the rather complex models that 
have been successfully treated with a blending syn- 
thesis method, as well as illustrating the choice of trial 
functions (planar solutions) and demonstrating the 
advantages of axially discontinuous single-channel 
synthesis. They considered the thermal reactor model 
shown in Fig. 2, represented by two-group diffusion 
theory. Nine radial trial functions were considered. 
Trial functions 1 to 5 were one-dimensional solutions 
representative of radial slices through the correspond- 
ingly numbered axial compositions. “Transition” trial 
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Fig, 2. PWR reactor model for two-dimensional synthesis test 
calculations.° 


functions 6 to 9 were one-dimensional radial solutions 
obtained using the averaged material constants of axial 
compositions 1—2, 2—3, 3—4, and 4—5S, respectively. 
Table 3 describes the trial functions used in two 
discontinuous synthesis solutions (1.a and 1.b), to- 
gether with the axial zone within which each set of 
trial functions was used, the motivation being to retain 
in any axial zone only those trial functions representa- 
tive of that zone and adjacent zones. 
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Table 3. Sets of Trial Functions for Discontinuous 
Synthesis Solutions” 








Test Fast trial Thermal trial Axial mesh 
problem functions functions range 
l.a 1,2 £2 0-9 
2.3 2,3 10-24 
3,4 3,4 25-38 
4,5 4,5 39-48 
1.b i52 1,2,6 0-9 
253 253,71 10-24 
3,4 3,4,8 25-38 
4,5 4,5,9 39—48 





Two continuous synthesis solutions were also 
performed, using a single set of trial functions over the 
entire axial domain. Continuous synthesis solution 1.c 
used trial functions | to 5 in both energy groups, 
resulting in 10 equations (two groups X five trial 
functions) to be solved, as compared to four equations 
(two groups X two trial functions) and five equations 
(two fast and three thermal trial functions) for 
discontinuous solutions 1.a and 1.b, respectively. Con- 
tinuous synthesis solution 1.d used trial functions 2 
and 3 in both energy groups, resulting in four 
equations to be solved. Galerkin weighting was used in 
all cases. 
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Thermal-flux distributions computed with the four 
synthesis calculations are compared with the results of 
a direct two-dimensional finite-difference solution in 
Figs. 3 to 6. Of particular interest is the superiority of 
synthesis solution 1.a relative to synthesis solution 1.d, 
both of which required the same core storage and 
computational effort to solve the synthesis equations, 
thus illustrating the advantage of using the discontinu- 
ous synthesis method. However, the computational 
effort needed to obtain trial functions was much 
greater for solution 1.a than for solution 1.d. The 
continuous synthesis solution 1.c is as good as or better 
than either of the discontinuous synthesis solutions. 
The former (1.c) required less effort in trial-function 
calculation and more effort in the solution of the 
synthesis equations (plus more core storage) than the 
latter. 

Depletion problems have been treated with syn- 
thesis by first doing depletion calculations on planar 
models corresponding to different axial zones. The 
planar flux solutions obtained at the different de- 
pletion steps were then used as trial functions in a 
three-dimensional synthesis depletion. 

Larsen?* recently examined the idea of recomput- 
ing trial functions at selected points during the 
depletion, thus eliminating the need for the planar 
depletion calculations. He reported several continuous 
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Fig. 3. Thermal-flux distribution: comparison of spatial synthesis with direct calculations.° 
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Fig.5 Thermal-flux distribution: comparison of spatial synthesis with direct calculations.” 


synthesis calculations for a two-dimensional two-group 
thermal reactor model representative of the Yankee 
pressurized-water reactor (PWR), as described in 
Table 4. His results are compared with the results of a 
direct two-dimensional solution (DBU) in Fig. 7. Re- 
computation of the trial functions resulted in a 
significant increase in accuracy, and there was no 
apparent advantage in using adjoint, rather than 
Galerkin, weighting. 


Although the majority of the applications of the 
blending synthesis method have been to thermal 
reactors, where the fine structure in the spatial flux 
distribution is quite detailed, some recent work?*:?° 
has demonstrated the utility of the method for 
fast-breeder reactor applications as well. Pilate et al.?* 
have made continuous synthesis calculations for three- 
dimensional models of the experimental facilities 
ZPR-3/48 and SNEAK-2C, using a five-group dif- 
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Fig. 6 Thermal-fiux distribution: comparison of spatial synthesis with direct calculations.’ 


Table 4 Continuous Synthesis (SYNTRON) 
Calculations for PWR Depletion? * 





Case Description 





{ Two trial functions in fast group and three 
in the thermal group; no trial-function recal- 
culation; Galerkin weighting 
il Same as case I but with trial-function recal- 
culation at 750 MWd/ton of uranium 
Il Same as case I but with trial-function recal- 
culation at 750 and 2000 MWd/ton of uranium 


IV Same as case III but with adjoint weighting 





fusion-theory representation, which have compared 
favorably with the results of a direct three-dimensional 
finite-difference calculation. The synthesis calculations, 
including the trial-function computations, were a 
factor of 10 or better quicker than the corresponding 
direct calculations. 

The geometrical model used for the SNEAK-2C 
calculation is shown in Fig. 8, where MASURCA and 
SNEAK refer to different fuel types and R2 is a radial 
blanket. Eight different planar shapes were considered 
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as trial functions, as indicated in Table 5. The first four 
corresponded to planar slices through the core and 
axial blanket regions with and without the control rod 
present, and the second four were auxiliary shapes, 
similar to Green’s functions (see the next subsection), 
which emphasized certain important regions. 

The Ker calculated for the reference rod-out and 
the rod fully inserted cases by the continuous synthesis 
calculations (KASY) are compared with the results of 
direct three-dimensional calculations (TRITON) and 
direct two-dimensional calculations through the core 
midplane with an estimated axial buckling (DIXY) in 
Table 6. The synthesis calculations were quite adequate 
for predicting Kegg and were remarkably accurate for 
predicting the worth of the control rod. Percent of 
errors in fission rate (absorption rate in the control 
rod) for the synthesis calculation, relative to the direct 
three-dimensional calculation, are plotted for z—y 
plane at y = 71.7 cm in Fig. 9. 

Luco?® compared continuous synthesis and direct 
finite-difference solutions for a two-dimensional model 
of a 1000-MW(e) LMFBR with middle-of-cycle (MOC) 
composition, represented by six-group diffusion 
theory. He used five different radial flux solutions as 
trial functions: trial function 1 was a radial solution 
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Table 5 Two-Dimensional Trial Functions?* for 
Synthesis Calculation of SNEAK-2C 




















Trial 
function Planar (x,y) configuration 
1 Plane at z = 57.35 cm with tantalum rod, B? = 1.16 x 10°*/cm? 
2 Plane at z = 57.35 cm without tantalum rod, B? = 1.16 x 10°3/cm? 
3 Plane at z = 20.00 cm with tantalum rod, B32 = 0 
4 Plane at z = 20.00 cm without tantalum rod, B? =0 
> Z, SNEAK fuel in rod position, tantalum absorber elsewhere, B3 = 0 
6 Z, SNEAK fuel in core zone 1, tantalum absorber elsewhere, B2 = 0 
7 Z, SNEAK fuel in core zone 2, tantalum absorber elsewhere, B? =0 
8 Blanket material in radial blanket, tantalum elsewhere, B =0 
Table 6 Criticality and Control-Rod-Worth Calculations 
for SNEAK-2C: Comparison of Spatial Synthesis 
(KASY) and Direct (TRITON) Solutions? * 
: Kett 
Methods Reference Rod in Ak (10°74) 
TRITON 0.9770 0.9600 170 
(DIXY)* (0.9813) (0.9639) (174) 
KASY+ 
One trial function (2) 0.9764 (1) 0.9594 170 
Two trial (2,4) 0.9768 (1,3) 0.9598 170 
functions (2,6) 0.9768 (1,6) 0.9598 170 
(2.7) 0.9767 
(2,8) 0.9766 
(1,5) 0.9595 
Three trial (2,4,6) 0.9769 (1,3,6) 0.9599 170 
functions (2,4,8) 0.9768 (1,3,8) 0.9598 170 





(1,2,3) 0.9598 





*Two-dimensional results obtained with an estimated B?. 
+The numbers of the trial function used are indicated in parentheses. 


through the core region; trial functions 2 to 4 were 
radial solutions through the three different axial 
blanket regions; and trial function S was a radial 
solution through the axial reflector, with an external 
volume source distributed as the fission source from 
the radial solution corresponding to the adjacent axial 
blanket region. 

Synthesis calculations were performed using the 5, 
3 (core, middle axial blanket, and reflector), and 1 
(core) trial functions, with Galerkin weighting. The 
results are compared to the results of a direct two- 
dimensional finite-difference calculation (DIF2D) in 
Table 7. (The designations C, AB, and RB indicate 
core, axial blanket, and radial blanket regions, re- 
spectively, and BR indicates the breeding ratio.) The 
synthesis calculations were quite accurate. 
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Luco performed another calculation using the same 
trial functions on the same model, but with a compo- 
sition corresponding to beginning-of-cycle (BOC), to 
test the possibility of using a single set of trial 
functions (calculated for one composition) for an 
entire fuel-cycle calculation. Results of the synthesis 
and a direct two-dimensional finite-difference calcula- 
tion are compared in Table 8. 


Green’s Functions Trial Solutions 


Green’s functions trial solutions were introduced 
by Dougherty and Shen?® and subsequently applied by 
others.2°-3! These authors suggested the use of flux 
solutions to problems in which the fission source was 
set to zero except in a limited region. By consecutively 
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Fig.9 Error (%) distribution for synthesized fission rate in SNEAK-2C. (Absorption rate in control rod.) 


Table 7 Spatial Synthesis Calculation for a Two-Dimensional LMFBR Model at Middle of Cycle? ° 





P ’ Pmax. 
— location, cm 
arbitrary 
Calc. Keff BR* Cl C2 ABI AB2 RBI RB units r z 





Regional 


Breeding Ratios 





DIF2D 1.003480 =1.37081 0.50449 = 0.28833 0.24182 =0.13552 0.11885 0.05181 
> LFF 1.003058 1.37429 0.50592 0.28708 0.27193 0.13376 0.14979 0.05583 
3 TF 1.003015 1.38121 0.50634 0.28673 0.24187 0.13375 0.14980 0.06272 
1 TF 1.002440 1.39521 0.50886 0.28509 0.23832 0.13604 0.14903 0.07796 





DIF2D 1.003480 


J EF 1.003015 
1 TF 1.002440 


1.37081 0.5071 0.3816 0.0515 0.0235 0.0311 0.0052 0.75254 
5 TF 1.003058 1.37429 0.5087 0.3798 0.0514 0.0233 0.0312 0.0055 0.75798 
1.38121 0.5092 0.3793 0.0514 0.0234 0.0312 0.0055 0.76197 
P.39S2E OSEES 0.3773 0.0499 0.0249 0.0310 0.0054 0.76933 


Regional Power Fractions 


1.266 3.810 
1.266 3.810 
1.266 3.810 
1.266 3.810 





*BR = breeding ratio. 
+TF = trial function. 


performing such calculations for all regions in a planar 
slice through a reactor, a set of trial functions was 
constructed, each of which was analogous to a Green’s 
function for the corresponding region. Trial functions 
5 to 8 in Table 5 are examples. Such trial functions 
should be well-suited for describing changes taking 
place separately in the different regions. This means of 
choosing trial functions has worked well for coupled- 


core reactor models*! but has been less successful for 
the analysis of more compact cores.” ? 


Multichannel Blending Methods of Spatial Synthesis 


Multichannel blending methods, as introduced by 
Wachspress and Becker® on a variational basis, have 
been extended to time-dependent problems by 
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Table 8 Spatial Synthesis Calculation for a Two-Dimensional LMFBR Model at Beginning of Cycle 
Using Middle-of-Cycle Trial Functions? ° 





Pmax Pmax. 
arbitrary _location, cm_ 
Calc. keff BR os | C2 ABI AB2 RBI RB2 units ar z 





Regional Breeding Ratios 


DIF2D 1.042968 1.46050 0.46276 0.28740 0.27912 0.17329 0.17703 0.08090 
Sar 1.040214 1.42670 0.47620 0.28065 0.28636 0.16808 0.14574 0.06966 





Regional Power Fractions 


DIF2D 1.042968 1.46050 0.4858 0.4112 0.0398 0.0238 0.0324 0.0069 
ey | 1.040214 1.42670 0.4998 0.4016 0.0407 0.0233 0.0286 0.0061 





0.69292 1.266 3.810 
0.69032 1.266 3.810 








Stacey®? and rederived on a weighted-residual basis by 
Fuller et al.2* This formalism allows the planar trial 
functions to be combined differently in different 
regions (channels) of the plane (or the use of different 
trial functions in different regions of the plane), which 
in essence allows the shape of the trial functions to be 
modified on a gross-region basis during the calculation. 
Such additional flexibility implies a greater accuracy 
and/or a greater range of applicable problems for a 
given set of planar solutions as trial functions, relative 
to the normal single-channel blending method, as has in 
fact been demonstrated by several tests.*?~3° 

However, numerical difficulties have been en- 
countered with the multichannel formalism which have 
hindered its implementation. The channel-to-channel 
coupling terms involve integrals over the mutual 
surface, and results have proved to be quite sensitive to 
numerical evaluation of these integrals. Moreover, the 
transport cross section is buried within these coupling 
terms following a series of matrix multiplications and 
inversions, making it burdensome to incorporate 
changes in the transport cross section (e.g., arising from 
density changes during a transient) into a calculation. 
Recent ideas on the use of overlapping polynomial 
domain functions to derive multichannel synthesis 
equations®” eliminate these two difficulties, as well as 
provide a formalism that can be treated with existing 
single-channel codes with only minor modifications in 
the processing of the trial functions. The multichannel 
blending method should reemerge as the most 
promising of the spatial synthesis methods. 


Synthesis of Reactivity Coefficients 


Synthesis of reactivity coefficients can be per- 
formed as an adjunct calculation, once the synthesized 
flux and adjoint distributions are constructed, by 
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performing the appropriate multidimensional perturba- 
tion integrals. Methods also have been developed for 
synthesizing three-dimensional reactivity coefficients 
directly?*'3° when it is not desirable or practical to 
construct explicitly the three-dimensional flux and 
adjoint distributions. 


Anomalies 


On occasion, anomalies have arisen in the appli- 
cation of spatial synthesis blending methods.*°’*? In 
addition, it has been possible to construct synthesis 
problems in which anomalous results were obtained for 
certain choices of cross sections and trial or weighting 
functions.*?-** These situations occur because the 
synthesis approximation, in general, fails to preserve 
certain positivity properties associated with the direct 
finite-difference multigroup diffusion equations. 
Although that eigenvalue of the latter which is associ- 
ated with the everywhere nonnegative (fundamental) 
eigensolution is always real and larger than the real part 
of any other eigenvalue, which ensures that iterative 
solution schemes that converge on the largest eigen- 
value will converge to the fundamental eigensolution, 
this is not necessarily true of the synthesis approxima- 
tion. Consequently it is possible for an iterative 
solution of the synthesis equations to converge to a 
solution that leads to a synthesized flux with negative 
components (harmonic) if the associated eigenvalue is 
larger than the eigenvalue corresponding to the funda- 
mental eigensolution or not to converge at all if the 
eigenvalue with largest real part is complex.*! 


This situation represents a serious but not fatal 
flaw in synthesis approximations. The danger is not so 
much in the truly anomalous case, where a negative 
flux serves as a warning that something is amiss, which 
can be recognized and rectified by changing trial or 
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weighting functions. The real difficulty is that, in the 
vicinity of that combination of composition and 
trial/weighting functions which gives rise to an obvious 
anomaly, the synthesized result is quite sensitive to 
small changes in any of these quantities; inaccurate but 
not implausible results may be obtained. Practitioners 
have developed empirical tests to guard against this 
difficulty, but these are in the realm of black art and 
are rarely written down. 

The situation is not so bleak as the preceding 
paragraph would indicate. Of the hundreds, perhaps 
thousands, of synthesis calculations that have been 
performed for realistic reactor models using reasonable 
trial and weighting functions, anomalies or associated 
difficulties have been found only a handful of times 
(the author knows of only three). Thus users and 
potential users of synthesis methods should be aware 
of but not distracted by the possibility of encountering 
such difficulties. 


Polynomial Expansions in Space 


Polynomial expansions in space were sug- 
gested*®*° several years ago as a means for arriving at 
more accurate but more highly coupled coarse-mesh 
difference equations. Considering the polynomials, 
which were nonzero only over certain (generally 
overlapping) regions, as trial functions, the methods are 
very similar, if not entirely equivalent, to spatial 
synthesis methods. 

A resurgence of interest has been kindled by the 
recent discovery (by reactor physicists/mathematicians) 
that the method has been highly developed by workers 
in structural mechanics, where it is known as the 
finite-element method.6®° The method is presently 
being investigated by several workers*'~*” for applica- 
tion to the neutron diffusion and transport equations. 

In addition to yielding more accurate coarse-mesh 
equations than are obtained by the usual finite- 
difference approximation, which may or may not be a 
practical advantage for realistic reactor models in 
which description of the geometric detail requires a 
rather fine difference mesh, the finite-element method 
is convenient for the treatment of complex geometries. 
Elements (within which polynomial expansions are 
made) may be chosen in any shape or size (providing 
one is willing to work out the algebra) and may vary 
within a given model. Although this is also true, in 
principle, for the derivation of the ordinary finite- 
difference approximation, it has been exploited to a 
much greater extent in the finite-element approxima- 
tion (e.g., Ohnishi’s work>?’>*) and has come to be 
regarded as a unique characteristic of the latter. 
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Stacey*” recently suggested a generalized nodal 
approximation in which spatial trial functions of the 
ordinary type are multiplied by polynomial domain 
functions, which may be considered as an extension of 
the finite-element method in which the polynomials 
are modulated by one or more spatial functions which 
account for details in the spatial distribution within an 
element. Henry*® recently proposed the use of a single 
spatial trial function to modulate the polynomial 
within an element. 


Coarse-Mesh Rebalancing 


Coarse-mesh rebalancing as a means of accelerating 
the convergence of flux iteration procedures by using 
the most recent flux iterate as a single trial function, 
which could be used to collapse to a coarser spatial 
mesh, was introduced by Wachspress*? within the 
context of variational synthesis. A theoretical founda- 
tion for the method was developed by Frohlich** *°° 
for the case of disjunctive, or nonoverlapping, coarse- 
mesh regions. He showed that the method retained the 
positivity properties of the fine mesh finite-difference 
equations for subdomain weighting and, under special 
circumstances, also for Galerkin weighting. 

The method has been incorporated in many dif- 
fusion-theory and Sy codes. Gains in convergence 
attainable with coarse-mesh rebalancing in the 
GAMBLE-5 diffusion theory code are shown in 
Table 9, as reported by Frohlich®' for a typical 
high-temperature gas-cooled reactor (HTGR) applica- 
tion. Impressive improvements were also reported for 
TRIGA applications. 

Nakamura®?~°* suggested the use of overlapping 
coarse mesh, obtained by superimposing overlapping 
domain functions on the most recent flux iterate, for 
rebalancing. Using “tent” domain functions, he demon- 
strated the superiority of the overlapping to the 


Table 9 Efficiency of Variational Rebalancing for a 
Typical HTGR Two-Dimensional Diffusion Theory 








Problem® ! 

Patch Execution 

Maximum Ket power time, min 

Rebal- _Itera- relative error, error, (UNIVAC 

ancing tions flux change % %o 1108) 

No 70 0.01 0.02 84 6.0 
No 140 0.003 0.0004 2.2 12.0 
Yes 70 0.005 0.0004 2.0 8.2 
Yes 140 0.0003 0.00001 0.06 14.7 
No 300 0.00003 Reference case 26.6 
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disjunctive coarse-mesh rebalancing for a _ one- 
dimensional diffusion theory model,°? as shown in 
Fig. 10. Here 6 is the Chebyshev parameter, K is the 
number of coarse-mesh regions, and F, is the error in 
the source. The rebalancing interrupted the normal 
Chebyshev accelerated iteration every n iterations. 
Nakamura®? also considered a three-group two- 
dimensional diffusion theory model representative of a 
boiling-water reactor (BWR). He compared a coarse- 
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Fig. 10 Decay of error source for overlapping (pyramid) and 
disjunctive coarse-mesh rebalancing methods for a 
one-dimensional diffusion theory problem. ° ? 


mesh rebalanced (every five iterations) and a normal 
Chebyshev accelerated iteration and found the former 
to be much superior (see Fig. 11). 


Summary 


Nonblending methods have been applied success- 
fully to simple geometric configurations and, with 
iterative improvement of the trial functions, to rela- 
tively more complex geometries. However, for more 
realistic models, blending methods offer a distinct 
advantage and have been developed to a much greater 
extent. The single-channel synthesis method has been 
applied to quite complex geometrical models with 
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Fig. 11 Decay of source error for variational rebalanced and 
ordinary Chebyshev iteration procedures for a two-dimensional 
BWR diffusion theory calculation.° ? 


rather impressive results, although the problem of 
obtaining trial functions for “driven” regions, such as 
blanket or reflectors, or for regions where the flux is 
highly nonseparable requires more work. Multichannel 
synthesis, which partially alleviates these difficulties, 
has encountered difficulties of a numerical type which, 
however, are eliminated by a new overlapping formula- 
tion. Thus it is expected that the development of an 
overlapping multichannel synthesis offers the potential 
for significant improvement in spatial synthesis blend- 
ing methods. 

The use of polynomial expansions in deriving 
coarse-mesh difference equations has been moderately 
successful in the past and currently is being investi- 
gated systematically. Coarse-mesh rebalancing, a special 
case of nonblending synthesis, has resulted in quite 
impressive gains when used to accelerate the iterative 
convergence of the finite-difference equations. Further 
development effort in both areas seems warranted on 
the basis of the results obtained to date. 


SPECTRAL SYNTHESIS 


Those methods in which precomputed spectral 
distribution functions (trial spectra) are blended to- 
gether to obtain a (generally detailed) solution for the 
neutron-energy spectra are referred to as spectral 
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synthesis methods. The names energy-modal approxi- 
mation, space—energy synthesis, and overlapping group 
have also been used to describe such methods. Calame 
and Federighi* were the first to develop and apply 
these methods. Subsequent authors have developed the 
theory,°*°® and the methods have been successfully 
applied to cross-section interpolation and space- and 
time-dependent calculations in thermal and fast reac- 
tors. 


Broad-Group Cross-Section Interpolation 


Broad-group cross-section interpolation is a natural 
application of spectral synthesis methods, and it is not 
surprising that it was one of the first such applications. 
Differential nuclear data were averaged over two or 
more different flux spectra (corresponding to different 
compositions) to form two or more broad-group 
cross-section sets. When the averaging fluxes were 
considered as trial spectra, a spectral synthesis calcula- 
tion yielded the appropriate combining coefficients for 
a given composition. These combining coefficients 
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could then be used to interpolate directly among the 
cross-section sets. 

Calame et al.°? combined “hard” (high absorption) 
and “soft” (low absorption) thermal-neutron trial 
spectra calculated for media consisting of 1/v absorbers 
to calculate thermal-group diffusion theory parameters. 
They employed different trial spectra for each differ- 
ent moderator temperature considered. A comparison 
of interpolated (SPG) thermal-group cross sections 
with thermal-group cross sections averaged over a 
composition-dependent spectrum calculated for the 
thermal-energy region by the multigroup code SOFO- 
CATE is given in Table 10 for a representative PWR 
composition. The interpolated cross sections are 
generally within 1% of the directly calculated values, 
yet the time required for synthesis solution and 
interpolation was approximately ,, of that for the 
SOFOCATE calculation. 
Ombrellaro et al.7°’7' considered the interpolation 
of fast-neutron broad-group cross sections. Using hard 
and soft trial spectra, Ombrellaro and Federighi’® 


Table 10 Comparison of Synthesized (SPG) and Directly Calculated (SOFOCATE) Thermal-Group 
Constants for a Representative PWR Composition® ° 


Hydrogen N = 0.03136 
Oxygen N = 0.01568 
Zirconium N = 0.02067 


Aluminum N = 0.002982 
735U N = 0.0001342 
738 U N = 0.00000825 














68°F 300° F a 500° F 
SOFOCATE SPG SOFOCATE SPG SOFOCATE SPG 
Diffusion Theory Parameters 
Sa, cm 0.07642 0.07551 0.06600 0.06519 0.05977 0.05903 
vd Ff, cm?! 0.13640 0.13460 0.11750 0.11588 0.10620 0.10465 
D,cm 0.34040 0.34330 0.36540 0.36720 0.38610 0.38690 
Microscopic-Absorption Cross Section (o,) 

Unit 1/y 0.7345 0.7258 0.6474 0.6397 0.5932 0.5863 
Hydrogen 0.2438 0.2409 0.2149 0.2123 0.1969 0.1946 
Oxygen 0.1469 x 1073 0.1451 x 10° 0.1295 x 1073 0.1279 x 10° 0.1186 x 10-3 0.1172 x 10-3 
Zirconium 0.1322 0.1306 0.1165 0.1151 0.1068 0.1055 
Carbon 0.2350 x 10°? 0.2322 x 10°? 0.2072 x 10°? 0.2047 x 10°? 0.1898 x 10°? 0.1867 x 107? 
Aluminum 0.1689 0.1669 0.1489 0.1471 0.1364 0.1348 
tl) 0.4882 x 10% 0.4824 x 10° 0.4202 x 10° 0.4150 x 10? 0.3797 x 10° 0.3750 x 10° 
eet) 0.2005 x 10! 0.1981 x 10! 0.1767 x 10! 0.1746 x 10! 0.1619 x 10! 0.1600 x 10’ 
Ree 0.2306 x 107 0.2227 x 10’ 0.2128 x 10’ 0.2071 x 10’ 0.1940 x 10’ 0.1895 x 10’ 
a7 %) 0.4316 x 10° 0.4271 x 10? 0.3796 x 10° 0.3763 x 10° 0.3482 x 10° 0.3456 x10? 
Hydrogen trans- 

port 0.2938 x 10? 0.2892 x 10? 0.2590 x 10? 0.2551 x 10? 0.2370 x 10? 0.2338 x 10? 

Microscopic Fission Production Cross Sections (var) 

aa0) 0.1016 x 10% 0.1003 x 104 0.8756 x 10° 0.8635 x 10° 0.7909 x 10° 0.7798 x 10° 
Ul 0.1004 x 104 0.9927 x 10% 0.8865 x 10% 0.8763 x 10° 0.8152 x 10° 0.8055 x 10° 
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interpolated fast broad-group cross sections (for repre- 
sentative PWR compositions) ‘which compared 
favorably with broad-group cross sections obtained by 
averaging over a multigroup spectrum calculated for 
the composition of interest. 

Ombrellaro and Snyder’! applied spectral syn- 
thesis to interpolate one- and three-group constants for 
a range of LMFBR core mixtures. They calculated 
50-group trial spectra for two extreme compositions 
(denoted by trial spectra 1 and 2 in Table 11), which 
were used to collapse to one- and three-group cross- 
section sets. Then a range of intermediate composi- 
tions, two of which are given in Table 11, was 


Table 11 Compositions for LMFBR Spectral Synthesis 








Calculations 7! 
Trial Trial 
Isotope Core 1 Core 2 spectra 1 spectra 2 
gn 0.012406 0.015531 0.009318 0.018640 
5©Fe 0.008962 0.010685 0.007293 0.012684 
oera) 0.004598 0.005905 0.003328 0.007214 
227 Pu 0.000861 0.000976 0.000732 0.001085 
PPE 0.002341 0.002791 0.001905 0.003313 
59 Ni 0.000922 0.001099 0.000750 0.001305 
241 Py 0.000067 0.000068 0.000051 0.000076 
240 Py 0.000344 0.000390 0.000293 0.000434 
23Na 0.012107 0.009880 0.014312 0.007556 
PFP1 0.000340 0.000426 0.000255 0.000511 








considered. Synthesis-extrapolated and directly calcu- 
lated (i.e., collapsed over a 50-group spectrum calcu- 
lated for that composition) group cross sections are 
compared in Table 12 for core 1. Similar results were 
obtained for core 2 and for cores in which the sodium 
was partially voided. 

These results indicate that the use of spectral 
synthesis to calculate composition-dependent combin- 
ing coefficients for the interpolation among broad- 
group cross-section sets is ari economic and relatively 
accurate means for obtaining composition-dependent 
broad-group cross sections. However, the ability of 
these methods to deal with the contribution of 
heavy-element resonance cross sections has not been 
adequately examined. 


Calculation of Spatially Dependent 
Thermal-Neutron Spectra by Spectral Synthesis 


Application of spectral synthesis to the calculation 
of spatially dependent thermal-neutron spectra was 
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Table 12 Comparison of Synthesized and Directly 
Calculated Few-Group Constants for a Typical LMFBR 
Composition’ ' 





Differences of synthesis values 
from exact values, % 











One-group Three-group 
hi: iti synthesis synthesis 
meters values i les ne P 

Groups 1—50 
D 0.20284+1 —0.05 0.003 0.02 0.21 
vir 0.61550-2 -0.12 -0.04 -0.02 0.21 
La 0.42335 -2 -0.11 -0.11 —0.04  -0.12 
Lt 0.56574 — 1 0.16 0.14 0.002 —0.025 
(p) 0.10764 + 3 0.12 0.04 0.02 -0.27 
VJ) 0.10917 + 2 0.06 0.05 0.03 —0.08 
WE rp) 0.66254+0 —0.003 0.0 -0.008 —0.06 

Groups 1-18 
D 0.23432 + 1 0.003 0.04 -0.01 0.13 
vir 0.65567 -2 —0.04 0.03 ~0.07 0.23 
La 0.32667 —2 —-0.03 0.004 —0.043 0.10 
Zt 0.53562 -1 —-0.01 —0.007 —-0.02 0.01 
{o) 0.71918+2 —-0.05 —0.01 0.08 0.05 
V) 0.84263+1  -0.04 0.03 0.07 0.18 


Groups 19-30 


D 0.14266 + 1 0.009 0.01 0.003 0.022 
vir 0.47841 —-2 —0.04 -0.04 -—0.04  -—0.05 
Xa 0.51509 —2 —0.07 —0.080  -—0.06 —0.10 
ry 0.58163 — 1 0.02 0.014 0.03 0.001 
(o) 0.33108 + 2 0.51 0.20 -—0.07 —0.83 
VJ) 0.23616 + 1 0.52 0.21 -0.07 —0.80 


Groups 31-50 


D 0.98474+0 —-1.2 —1.1 0.05 0.06 
vir 0.12465 -—1 —-0.76 -0.61 0.89 0.91 
La 0.19206 -1 -0.70 —0.56 0.77 0.79 
zy 0.11928 + 0 32, 29 0.52 0.49 
(p) 0.26150 + 1 0.39 —0.42 —0.72 —1.8 
V) 0.12876+0 -1.6 —1.5 —0.66 —1.7 





first carried out by Calame and Federighi.* They 
argued that the thermal-neutron spectrum in a two- 
region model consisted of two components: (1) neu- 
trons which became thermal in region 1, with a thermal 
spectrum characteristic of that medium, and thereafter 
diffused throughout space while retaining that spec- 
trum and (2)neutrons which became thermal in 
region 2 and thereafter diffused throughout space 
while retaining that characteristic spectrum. Thus they 
were led to blend infinite media spectra characteristic 
of the two regions, with the knowledge that one or the 
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other was appropriate far away from the interface 
between the two media and with the hope that a 
combination of the two infinite media spectra would 
be a good approximation in the vicinity of the 
interface. Calame and Federighi developed coupled 
equations for the spatially dependent combining co- 
efficients and successfully applied the method to the 
analysis of experiments. The method subsequently has 
been incorporated into a number of codes and used 
extensively, yielding results comparable in accuracy to 
those obtained by a more expensive multigroup treat- 
ment of the thermal-energy spectrum in a spatially 
dependent calculation. 

Extensions and refinements of this method, which 
is generally known as the two-overlapping-group 
method, have been suggested by several authors. 
Leslie’? considered the terms that coupled the two 
equations for the combining coefficients as “removal” 
cross sections from one spectral distribution to another 
and established a physical basis for their definition and 
computation. Breen and Yasinsky’* proposed a par- 
ticular choice of weighting spectra which partially 
decoupled the two equations for the combining co- 
efficients so that they could be solved with con- 
ventional two-group iteration schemes. Natelson’* 
employed the method to develop a two-overlapping 
thermal-group approximation for the discrete-ordinates 
representation of the transport equation, used a tech- 
nique similar to that of Breen and Yasinsky to simplify 
the equations, and obtained results that compared 
rather well with Monte Carlo calculations. Selengut’* 
examined the choice of weighting spectra, which 
traditionally are chosen as F (the energy variable) and 
unity, and suggested that \/E was a better choice than 
E. 


Fast Reactor Applications of Spectral Synthesis 


Fast reactor applications of spectral synthesis were 
first investigated by Stacey®® and by Storrer and 
Chaumont.’® Stacey considered a two-region model 
(core and blanket) representing a fast reactor critical 
facility. As trial spectra he chose: (1) a zero leakage 
spectrum characteristic of the core medium and (2) the 
spectrum resulting from solving a source problem in an 
infinite medium of the blanket material, using as the 
source the 73° U fission spectrum characteristic of the 
core medium. He demonstrated that the first trial 
spectrum was appropriate at the center of the core, 
that the second trial spectrum was a good approxima- 
tion deep in the blanket, and that a mixture of the two 
provided a reasonable approximation in the vicinity of 
the interface. 
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These early studies have been extended by several 
authors’ ’~®* in an attempt to establish the feasibility 
of spectral synthesis for fast reactor analysis, where the 
possibility of obtaining the necessary spectral detail by 
blending a few trial spectra promises significant 
economy relative to a direct calculation with the large 
number of groups which would be required for an 
equivalent multigroup spectral representation. Differ- 
ent choices of trial and weighting spectra have been 
compared, and the advantages of using different sets of 
trial or weighting spectra in different spatial regions has 
been examined. Neuhold and Ott’? and Neuhold®°® 
introduced the idea of weighting with reaction rates, 
which they found to be generally superior to Galerkin 
weighting and comparable to adjoint weighting. Errors 
in the synthesis calculation of the absorption rate in a 
slab model of an LMFBR are shown in Fig. 12. Because 
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Fig, 12 Error in the synthesized absorption rate in a 
one-dimensional LMFBR model for different weighting 
functions.* ° 


the purpose of this study was to assess different 
weighting spectra, the trial spectra were taken as the 
exact solutions at the center of the core and at the 
center of the blanket. The two weighting functions in 
the case designated “reaction rate” were the product of 
the core macroscopic fission cross section with the two 
trial spectra. For the case designated “‘multiple reaction 
rate,’ the two weighting functions used in the core 
were the products of the core macroscopic-absorption 
cross section with the two trial spectra, and the two 
weighting spectra used in the blanket were the prod- 
ucts of the blanket macroscopic-absorption cross 
section with the two trial spectra. “Galerkin” refers to 
the use of the trial spectra also as weighting spectra, 
and “adjoint” refers to the use of the exact adjoint 
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spectra at the center of the core and at the center of 
the blanket as weighting spectra. Best results with 
reaction-rate weighting were obtained by using differ- 
ent reaction rates when calculating different integral 
quantities, which becomes impractical. 

Several authors’?~°? have found that the use of 
different sets of trial and/or weighting spectra ap- 
propriate to the different regions of a reactor generally 
results in a better approximation than can be obtained 
by using the same number of trial spectra throughout 
the reactor. Neuhold and Ott’? and Kiguchi et al.*' 
found that, in regions where the spectrum was largely 
influenced by another region (e.g., a blanket), a 
spectrum calculated for the “driven” medium with a 
leakage source or a buckling representing the neighbor- 
ing “driver” region was a better trial spectrum than an 
isolated fundamental-mode solution for the driven 
medium. Neuhold and Ott’? demonstrated the im- 
portance of using an adequate fundamental-mode 
leakage representation for both core and blanket 
calculations. 

Kiguchi et al.°' found that the sodium-void worth 
in an LMFBR model was predicted more accurately by 
a spectral synthesis calculation than by a conventional 
few-group calculation requiring equivalent computing 
time. Stacey®* concluded that the spectral synthesis 
method was superior to few-group methods (requiring 
equivalent computing time) for treating the reactivity 
worth of sodium voiding and the interacting Doppler 
worth during a transient in a space-independent 
LMFBR model. Both authors used trial spectra calcu- 
lated for core regions with extremes in sodium con- 
centration to synthesize the spectrum for intermediate 
sodium concentrations. Stacey®* also demonstrated 
the ability of spectral synthesis to treat local sodium 
voiding in spatially dependent kinetics calculations. 


Multidimensional Geometries in Spectral 
Synthesis Equations 


Multidimensional geometries pose special problems 
for the solution of the spectral synthesis equations. 
Although the spectral synthesis and multigroup equa- 
tions are formally identical, the coupling among 
equations is much greater for the former, with the 
result that the diffusion coefficient and removal 
matrices, which are diagonal for the multigroup equa- 
tions, and the scattering matrix, which generally is 
lower triangular for the multigroup equations, are all 
full for the spectral synthesis equations. Thus the 
iterative solution techniques, which have been de- 
veloped for the multigroup equations and which take 
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advantage of the special nature of these equations, 
cannot be used directly for the solution of the spectral 
synthesis equations. This is particularly true for the 
inner iterations, which are used to invert matrices, and 
the acceleration of the outer iterations. 

Although there are two possible means of resolving 
this difficulty, modification of the multigroup iteration 
schemes together with recasting of the spectral syn- 
thesis equations or the development of new iteration 
strategies especially for the spectral synthesis equa- 
tions, only the former alternative has been explored. 
Robkin®® and Lorenzini and Robinson®’ suggested a 
series of transformations which reduced the spectral 
synthesis equations to a form that can be solved by 
multigroup codes. Such transformations basically in- 
volved the diagonalization of the diffusion-coefficient 
matrix, followed by the addition to the fission source 
term of those parts of the scattering and removal 
matrix which do not conform to the multigroup 
structure. Lorenzini and Robinson®” were able to solve 
two-dimensional LMFBR spectral synthesis problems 
with a standard multigroup code, although they gen- 
erally found that more iterations were required than 
for multigroup problems with the same number of 
unknowns. Numerical instabilities associated with the 
complete coupling in the scattering matrix (they used a 
code that allowed upscattering and hence a full 
scattering matrix) were encountered, and modifications 
of the diagonal terms were introduced to suppress this 
difficulty. Although they compared spectral synthesis 
and few-group approximations, their results were in- 
conclusive as to which was superior. They introduced 
the idea of using one-dimensional calculations to 
generate trial spectra appropriate to different regions in 
a two-dimensional model. 

Greenspan®*® also applied spectral synthesis to 
two-dimensional LMFBR models, using a formulation 
similar to that of Lorenzini and Robinson to enable the 
equations to be solved with a standard multigroup 
code. His model consisted of two core regions, axial 
and radial blankets, and a radial reflector. Trial spectra 
(24 group) were: (1) a fundamental-mode solution for 
the dominant core material, (2) a source solution for 
the axial blanket material utilizing the leakage 
spectrum from the core trial spectrum as a source, and 
(3) a source solution for the reflector material utilizing 
the leakage spectrum from the blanket trial spectrum 
as the source. Greenspan employed subdomain (in 
energy) weighting whereby the 24 groups were dis- 
tributed among three energy intervals and each 
weighting spectrum was unity within one energy 
interval and zero within the others. For comparison, he 
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performed several few-group calculations in which the 
few-group constants were collapsed over the spectrum 
of the dominant core region (i.e., the first trial 
spectrum). He found that the spectral synthesis calcula- 
tions required on the order of 10% greater computing 
time than the few-group calculations with the same 
number of unknowns (i.e., in which the number of trial 
spectra and few groups were equal) but that consid- 
erably more few groups than trial spectra were required 
to obtain comparable accuracy, so that the spectral 
synthesis approximation was superior to the few-group 
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approximation in this application. Computing time 
requirements and important reactor parameters are 
compared with the results of direct 24-group calcula- 
tions in Tables 13 and 14. 


Iterative Improvement of Trial Spectra 


Iterative improvement of trial spectra has been 
investigated by several authors?®°? Two general 
methods have been considered; successive iteration 


between spectral and spatial synthesis®* "°° *?° and the 


use of approximate solutions together with a single- 
group diffusion kernel to arrive at zero-dimensional 
equations from which an approximation to the 
spectrum at a point in a reactor can be obtained.?' °°? 
Toivanen,®’ considering a small intermediate- 
spectrum reactor, found that asymptotic spectra repre- 
Time, sentative of the different regions in the reactor were 

k ett sec inadequate trial spectra. He found that iteration 
between spatial and spectral synthesis, as suggested by 
Stacey,°* significantly improved the trial spectra and 


Table 13 Comparison of k, ¢¢ and 
Computing Times for a Two- 
Dimensional LMFBR Problem® ® 








Reference Case 


24 group 1.05883 414 hence the accuracy of the calculation. A similar idea 
3 TF synthesis 1.05911 109 ; cued ty Lencaibald”™ ie oie inehea’ 
Seve 1.07787 84 was pursue 7 ancefie in app ying spectra 
4 group 1.07064 102 synthesis to a fast reactor model. Lancefield also 
8 group 1.06412 149 suggested the use of homogeneous functionals of the 


flux, such as reflector reaction rates, in deriving the 


Sodium Voided in Dominant Core synthesis equations so as to improve the accuracy of 


24 group 1.05447 430 the flux solution in the reflector. This latter idea is 
3 TF synthesis 1.04706 99 similar to Neuhold’s®® use of multiple reaction-rate 
3 group 1.07446 88 weisbtins 
4 group 1.06665 105 eae 





Cockayne and Ott,?' based on studies of a 
two-region (core-blanket) LMFBR model, concluded 
that the actual spectra occurring at selected points in a 


*Same trial functions and collapsing 
spectrum used as in the reference case, 


Table 14 Power Fractions and Breeding Ratios for a Two-Dimensional 








LMFBR Problem®® 
Axial Axial Radial 
Core | Core 2 blanket 1 blanket 2 blanket Total 
Power Fraction 
24 group 0.6572 0.3233 0.0080 0.0032 0.0083 
3 TF synthesis 0.6609 0.3162 0.0105 0.0040 0.0084 
3 group 0.6458 0.3239 0.0130 0.0048 0.0125 
4 group 0.6513 0.3776 0.0089 0.0035 0.0087 
8 group 0.6541 0.3227 0.0096 0.0037 0.0099 
Breeding Ratio 
24 group 0.5080 0.1496 0.2124 0.0691 0.1932 1.1323 
3 TF synthesis 0.5096 0.1479 0.2178 0.0706 0.1908 1.1367 
3 group 0.5032 0.1547 0.1877 0.0626 0.1805 0.0886 
4 group 0.5048 0.1535 0.1930 0.0642 0.1828 1.0983 
8 group 0.5077 0.1512 0.2035 0.0647 0.1891 0.1162 
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reactor made the best trial spectra. They examined 
several methods of approximating such spectra. The 
fundamental-mode spectrum of the core material was 
found to be a very good approximation to the actual 
spectrum over much of the core, but the fundamental- 
mode spectrum for the blanket material was found to 
be a very poor approximation to the actual spectrum in 
the blanket. The blanket spectrum, calculated using as 
a source the fundamental-mode core leakage spectrum, 
was found to be a good approximation to the actual 
spectrum over a portion of the blanket and resulted in 
a relatively accurate synthesis solution when used as a 
trial spectrum. These authors developed an approxi- 
mate procedure for computing the spectrum at a point 
in a reactor model which consists of three steps: 
(1) calculation of an approximate flux from a two- 
spectra synthesis problem using two readily available 
spectra; (2) calculation of special “group constants” at 
the specific spatial location using as weight functions 
the approximate flux from step 1 and a single-group 
diffusion kernel; and (3) solution of the resulting 
zero-dimensional equation for the spectrum at the 
point, which is then used as a trial spectrum. Use of 
this procedure in one LMFBR problem significantly 
reduced the error in blanket absorption rates (from 
1.27 to 0.6%) and fission rates (from 3.64 to 0.17%) 
relative to the results of step 1. This procedure may be 
limited, in practice, to one-dimensional problems be- 
cause of complexities associated with the diffusion 
kernel in multidimensional geometries. 


Summary 


Spectral synthesis has been demonstrated to be an 
economical and accurate approximation to a detailed 
multigroup representation of the neutron spectrum for 
many applications and generally to be superior to the 
conventional few-group reduction of the multigroup 
representation. If detailed multigroup trial spectra are 
blended, the spectral synthesis method has the po- 
tential for representing the spectrum in great detail 
with a relatively simple calculational model. As was the 
case with spatial synthesis, the spectral synthesis 
approximation does not retain the positivity properties 
of the multigroup equations, and anomalous results can 
be obtained. 

The use of spectral synthesis to interpolate among 
broad-group cross-section sets averaged over different 
fine-group spectra has been successfully demonstrated 
for the thermal- and fast-neutron energy regimes. 
However, treatment of the composition and 
temperature dependence of heavy-element (e.g., 
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uranium) narrow-resonance cross sections should be 
developed and evaluated to a greater degree. 

Spatially dependent spectral synthesis (overlapping 
thermal-group) methods have been successfully applied 
to the thermal-neutron regime for many years. Applica- 
tions to the fast-neutron-energy regime and to fast 
reactors have been more _ recent. However, the 
adequacy and economy of spectral synthesis calcula- 
tions for fast reactors seem to be relatively well 
established. The fact that the spectral synthesis equa- 
tions (suitably modified) can be solved by standard 
multigroup codes enhances the attractiveness of the 
method. The calculation of trial spectra representative 
of “transition” regions near region interfaces substan- 
tially improves the accuracy in these transition regions. 
Two areas in which more work is needed are: 
(1) adaptation of standard multigroup iteration 
schemes and of the spectral synthesis equations to 
achieve an efficient means for solving the latter in 
multidimensional geometries and (2) development of 
economical means for computing trial spectra charac- 
teristic of “transition” regions and other regions where 
the asymptotic spectrum characteristic of the medium 
is a poor approximation to the actual spectrum. 


TRANSPORT-THEORY APPLICATIONS 
OF SYNTHESIS METHODS 


In the preceding sections the use of synthesis 
methods to solve the multigroup diffusion or P; 
equations has been emphasized. In this section, three 
types of applications of synthesis methods to neutron- 
transport theory are reviewed. The first such applica- 
tion is the derivation of approximations to the 
transport equation in which the continuous angular de- 
pendence is discretized. The second application is to 
the solution of these discretized equations, using 
methods similar to those which have been discussed in 
previous sections. Finally the application of spatial 
synthesis methods directly to the transport equation, 
preceding the discretization of the angular dependence, 
is discussed. 


Discretization of Angular Dependence 


Discretization of the angular dependence is the 
normal first step in seeking a solution to the neutron- 
transport equation. Expansion in spherical harmonics 
and direct discretization, leading to Py and discrete- 
ordinates approximations, respectively, are traditional 
methods in reactor physics. Selengut® showed that, 
when the first two Legendre polynomials were con- 
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sidered as angular trial functions in a variational 
derivation, the resulting synthesis equations were the 
P, equations. This insight precipitated the develop- 
ment of a number of discrete approximations which 
were variants on the Py or discrete-ordinates approxi- 
mations. 


Pomraning and Clark?? derived an improved form 
of diffusion theory, which had an asymptotic trans- 
port-theory diffusion coefficient, by making a two- 
term expansion of the angular dependence of the 
neutron flux. Numerical comparisons based on a 
thermal-group cell problem indicated that his improved 
diffusion theory was comparable in accuracy to P3 
theory. Palveri-Fontana and Amster?* also derived an 
alternate diffusion theory, featuring a different defini- 
tion of the diffusion coefficient, based on a double Po 
expansion of the angular distribution and utilizing a 
discontinuous trial-function variational formulation. 


Kaplan et al.?*° used circles and ellipses as angular 


trial functions to reduce the transport equation to an 
angular discrete approximation. They noted that dis- 
crete-ordinate approximations could be derived in this 
manner by choosing the angular trial functions as 
constant within each solid-angle element. Schreiner and 
Selengut?® proposed a similar derivation of discrete- 
ordinate equations, using linear polynomials in the 
cosine of the angular variable within each solid-angle 
element. 


Natelson?” suggested a strategy for deriving a 
discrete approximation by using at each spatial point 
the angular distribution (in a high-order discrete- 
ordinate approximation) obtained from a reference 
calculation as an angular trial function. He also 
suggested some approximate strategies for obtaining 
the angular trial function when a high-order discrete- 
ordinate approximation was impractical. Using these 
angular trial functions (one at each spatial mesh point 
but different trial functions at different mesh points) 
obtained for a reference or “representative” problem, 
he derived diffusion-theory-type approximations that 
were then applied to similar problems. He obtained 
impressive results for this space-angle-synthesis (SAS) 
approximation for a variety of realistic problems. One 
such problem is depicted in Fig. 13, and the results 
obtained with SAS and P, diffusion theory are 
compared to a high-order discrete-ordinate (D,,) 
solution in Table 15. 


Several authors?*!°° have examined the varia- 


tional synthesis derivation of discrete-ordinate equa- 
tions. Jauho and Kalli?® proposed two methods for 
choosing the set of directions and weights to be used in 


Weston M. Stacey, Jr. 229 












































j 
I (Fuel) 
Eg. = 0.3169 cm’! 
I =_ 
Z_ = 0.6033 cm! 3 
bs ee E 
o ro) ° 
S, = 0.015 so) 5 rs 
Ss] § g 
s) bc 
) “ 5 
> 1 
DT = 10cm NO 
Sy = 0 
ble | 
0 (Fuel) 
Zs, 7 0.3694 cm’! E 
§ 
D,_ = 0.5693 cm’! = 
Oo o 
Sy = 0.015 
i a ais. cl st Sle ashi a ig SC 3 at- 
a 
a15 
tL} jsy° 
| , 
I (Structure and water) a 
Ze = 0.5314 cm’! VII (STRUCTURE) Ze 
7 z = 0.1726 cm"! > 2 
Z,__ = 0.5443 cm"! “Syn uw & 
~ E_ = 0.1804 cm! Ss 
Spy = 0.02 Tym = 
Syn = 0 z= 
Ww 
= 
oO 
re) 
IV (Water) = 
(5) 
Lg = Zg_ = 0.8729 cm! - 
Ly Vv Ww £ 
E,_ = 2,_ = 0.8908 cm’ ba 
Iv Vv N 
Sry = Sy = 0.06 | 
a 2.278 cm | | as 
(38 MESH INTERVALS) 
0.244 cm — 
0.146 cm +— 


Fig. 13 Heterogeneous PWR cell model for space-angle- 
synthesis calculations.° ” 


the standard Sy equations. Kaplan®® used angular trial 


functions that were constant within each solid-angle 
element, but oriented the solid-angle elements dif- 
ferently at adjacent spatial mesh points, in seeking 
discrete-ordinate equations that would be free of the 
troublesome “ray” effect. 

Natelson' °° found that a variational derivation of 
the standard discrete-ordinate approximation in two 
dimensions also yielded a relation between the weights 
and the directions but that use of this relation did not 
produce results significantly different from those of 
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the standard discrete-ordinate scheme. He also derived 
a modified discrete-ordinate approximation that con- 
sisted of pairs of coupled diffusion-type equations, 
with each increase in the order of the approximation 
adding another pair of equations. The potential ad- 
vantage of the modified approximation is its inherent 
ability to mitigate ray effects, thus to serve as a 
standard. 


Spatial Synthesis in Transport Theory 


Spatial synthesis in transport theory has developed 
along two lines. One procedure involves the blending of 
discretized transport-theory solutions of lower dimen- 
sional (in space) problems to obtain higher dimensional 
solutions.'°'"'°? A second procedure consists in 
reducing the spatial dimensionality by blending spatial 
trial functions, then carrying out the angular discretiza- 
tion on the spatially reduced equations.'°*! °° 

The first procedure is typified by the work of 
Natelson and Gelbard,'°? who blended two-dimen- 
sional discrete-ordinate solutions with combining co- 
efficients which depended on the third spatial variable 
but not on the angular variable. The equations that 
resulted for the combining coefficients were of the 
same form as in the case of diffusion theory synthesis, 
which facilitated implementation of the method in an 
existing diffusion-theory spatial synthesis code. This 
method has the potential for synthesizing angular 
distributions in the plane of the trial-function calcula- 
tion but has no provision for treating angular effects 
not in that plane (i.e., it can handle azimuthal but not 
polar angular effects). 

Cobb’s'°® work represents the most extensive 
development of the second procedure. He blended 
one-dimensional spatial shapes, with combining co- 
efficients which depended on the other spatial variable 
and the angular variable, to reduce a two-dimensional 
transport equation (one group) to a coupled set of 
one-dimensional transport equations. A _ discrete- 
ordinate approximation was then made for the latter. 
He used the multichannel synthesis approximation to 
allow different combining coefficients to be used in 
different spatial regions. A comparison with a direct S4 
two-dimensional calculation is given in Table 16 for a 
small 10-cm square fuel assembly with a 4-cm-square 
absorber in one corner. The two trial functions were 
one-dimensional Sq calculations along traverses that 
did and did not pass through the absorber. These 
results indicate that this type of synthesis approxima- 
tion can yield rather accurate results, even for this 
highly nonseparable problem, together with a consider- 
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able savings in computing time relative to the direct 
solution. 

Lancefield’s’?® work with spectral synthesis also 
utilized the second procedure. He used trial spectra to 
discretize the energy representation of the transport 
equation and then used a discrete-ordinates approxima- 
tion to solve the resulting equations. 


Summary 


The use of variational synthesis methods to dis- 
cretize the angular dependence of the neutron- 
transport equation not only provides a consistent 
means of deriving the standard (e.g., Sy, Py) approxi- 
mations but also provides a systematic procedure for 
deriving alternate approximations. Exploitation of this 
technique has been limited, but the results have been 
promising, suggesting that further effort in this area is 
warranted. 

Standard spatial synthesis methods, involving the 
blending of two-dimensional transport solutions with 
angle-independent combining coefficients, are in- 
herently incapable of treating angular effects not in the 
plane of the two-dimensional solutions. An extension 
of these methods to use angle-dependent combining 
coefficients would seem to be worthwhile. 

A second type of spatial synthesis, in which a 
spatial expansion is used to reduce the dimensionality 
of the problem prior to the discretization of the 
angular dependence, also has been employed with 
encouraging results. Further investigation is required to 
determine if this type of spatial synthesis has any 
advantage or disadvantage with respect to the standard 
type. 


VARIATIONAL PRINCIPLES 


Significant advances in flux synthesis methods have 
proceeded, almost without exception, from advances in 
the underlying variational theory. Although the entire 
theory of flux synthesis methods can be constructed 
with the method of weighted residuals, and all the 
advances can be derived therefrom without recourse to 
variational methods, the fact remains that variational 
theory has played a central role in the development of 
flux synthesis methods. Thus this section is devoted to 
a review of those developments in variational theory 
which are relevant to flux synthesis methods. Construc- 
tion of variational principles, extension of variational 
principles to admit a wider class of trial functions, and 
the transformations that relate different variational 
principles will be discussed. 
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Construction of Variational Principles 
for Neutron Transport and Diffusion Theory 


Construction of variational principles for neutron 
transport and diffusion theory has been, one suspects, 
largely a matter of inspiration. However, some effort 
has been devoted to the development of procedures for 
systematically deriving variational principles. 

Pomraning'®”"'!° and Lewins''':''? have sug- 
gested a method to derive variational principles for the 
estimation of some integral parameter. They adjoined 
the neutron transport/diffusion equations to the func- 
tional which defined a parameter of physical interest 
by means of Lagrange multipliers, in the way that 
constraints are treated in the formulation of a problem 
in the classical calculus of variations. It turned out that 
the Lagrange multipliers must satisfy the adjoint 
equations in order for the resulting functional to be a 
variational principle. 

Becker''? and Mikhlin''* suggested “‘least- 
squares” variational principles formed by “squaring” 
the neutron transport/diffusion equations and inte- 
grating over the independent variables. Such principles, 
which have the advantage of not requiring the use of an 
adjoint but the disadvantage of higher order dif- 
ferential operators in the final synthesis equations, 
have been little used. 

Selengut''® and Pomraning'®® developed varia- 
tional principles by considering a functional Taylor’s 
series expansion of the functional which defined the 
integral parameter of interest. Termination of the series 
with the quadratic term, plus the requirement of 
accuracy to second order in the trial function, led to 
the appropriate variational principle. 

Pomraning'®® related this latter derivation to 
perturbation theory and indicated that retention of 
higher order terms in the series led to higher order 
variational principles (i.e., principles for which errors in 
the trial function led to nth order (n > 2) errors in the 
value of the functional) of the type proposed by 
Kostin and Brooks.''® Devooght''? also derived 
higher order variational principles. These higher order 
principles do not seem to have been employed yet in 
the derivation of flux synthesis approximations. 


Extended Variational Principles 
in Flux Synthesis Methods 


Extended variational principles have played a major 
role in the development of more powerful and more 
flexible flux synthesis methods. A variational principle 
for the neutron transport/diffusion problem must not 
only have the appropriate balance equations as Euler 
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equations but must also embody the associated 
boundary, initial, and continuity conditions, either 
directly or indirectly through restrictions on the 
admissible class of trial functions. Significant advances 
have been made in embodying these associated “side” 
conditions directly in the variational principle, thus 
increasing the class of admissible trial functions and 
thereby the flexibility and power of the flux synthesis 
method. 

Several authors have considered the 
problem of embodying external surface boundary 
conditions directly in the variational principle, thus 
enabling the derivation of a synthesis approximation in 
which the trial functions were not constrained to 
satisfy the external boundary conditions. Such varia- 
tional principles were also employed in the derivation 
of appropriate external boundary conditions for the 
Py equations. 

Lewins'*® and Becker'*? noted that the varia- 
tional principle of Kaplan et al.'? rigorously admitted 
only flux trial solutions satisfying a known final 
condition and introduced an initial value term to 
replace this physically unrealistic condition by initial 
value conditions on the trial solutions. Pomraning' ?° 
and Yasinsky'® suggested initial and final value terms, 
which embodied flux initial conditions and adjoint 
final conditions directly in the principle, and temporal 
continuity terms, which embodied flux and adjoint 
temporal continuity directly in the principle. The 
former enabled physically realistic initial and final 
conditions to be derived for the flux and adjoint, 
respectively, combining coefficients in a synthesis 
approximation, and the latter enabled different sets of 
trial functions to be blended within different time 
intervals. 

Wachspress and Becker® incorporated flux and 
current spatial continuity conditions directly into a 
variational principle, thereby enabling the use of trial 
solutions that need not satisfy the physical conditions 
of flux and current continuity. Multichannel blending 
synthesis® and discontinuous synthesis? were direct 
results of this added flexibility. Initially it was held 
necessary to work with a variational principle for the 
P, equations in order to embody spatial continuity 
conditions. However, Pomraning'?® and Buslik'?! 
introduced appropriate variational principles for the 
neutron diffusion equation, and Lambropoulos and 
Luco!3? subsequently provided mathematical argu- 
ments as to the validity of such principles. Goldstein 
and Lancefield'*? have suggested a potentially useful 
theory in which the locations of the spatial interfaces 
are determined variationally. 


6,118-127 
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Transformations and Complementary 
Variational Principles in Transport Theory 


Neutron-transport theory can be characterized by 
many different equations that are related to each other 
by various transformations. To each characterization 
corresponds a different variational principle. It was 
recently shown that many of these variational 
principles are related among themselves by transforma- 
tions of two special types, canonical and involutory. 
Moreover, it has been demonstrated that the pair of 
extremum (minimum and maximum) variational 
principles that are related through an involutory 
transformation are complementary, i.e., the functions 
that minimize one maximize the other. 

Kaplan'?* drew an analogy between reactor 
physics and classical mechanics to demonstrate that the 
diffusion-theory variational principles of Selengut? and 
Kaplan> were related to the P; theory variational 
principle of Wachspress and Becker® in the same 
manner that Hamilton’s principle is related to the 
canonical integral. Kaplan and Davis'?* demonstrated 
the relation among four variational principles for the 
neutron-transport problem. Beginning with a gen- 
eralized form of one of Vladimirov’s principles,'*° 
which has the even-parity second-order form of the 
transport equation as its Euler equation, they applied a 
canonical transformation to obtain essentially Davis’ 
principle'?? F,, but in a form that did not contain 
adjoint functions. Then they applied an involutory 
transformation that resulted in a principle related to 
Pomraning and Clark’s principle’ F[¢], but with 
surface terms, which has the odd-parity second-order 
form of the transport equation as an Euler equation. 
Applying a canonical transformation to this involutory 
principle produced essentially Davis’ principle'** Fy, 
which, with F,, has the transport equation in coupled 
first-order form as Euler equations. Adding a variable 
surface term to Vladimirov’s principle resulted in a 
family of principles, each of which has canonical and 
involutory forms. One of the canonical forms turned 
out to be essentially Selengut’s principle.''® 
Kaplan'?7 extended the theory for canonical and 
involutory transformations to variational problems 
with higher derivatives. Pomraning'?® further ex- 
tended the theory to allow an arbitrary number of 
linear operators, including derivatives of all order, in 


the variational problem. 


Kaplan and Davis'** pointed out that the gen- 


eralized form of Vladimirov’s principle and the 
principle obtained after a canonical and an involutory 
transformation are complementary (for monoenergetic 
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transport theory) in the sense that the minimum value 
of the former equals the maximum value of the latter. 
Their suggestion that this fact might lead to a means 
for obtaining error bounds or figures of merit for 
approximate solutions was investigated by Yasinsky 
and Kaplan.'*? These latter authors applied the idea 
to obtain a figure of merit for synthesis approxima- 
tions in neutron-diffusion theory with some success 
but concluded that the computational effort involved 
rendered the procedure unfeasible for large multidi- 
mensional problems. 

Pomraning' *8*'*° used involutary transforma- 
tions to construct complementary variational principles 
for self-adjoint neutron-transport and diffusion prob- 
lems. He invoked these complementary principles to 
compute bounds for several quantities of interest in 
neutron-transport and reactor theory. 

Davis'*' formulated complementary variational 
principles for the monoenergetic one-dimensional 
neutron-transport equation. He utilized these com- 
plementary principles to show that even-order Py 
approximations converge monotonically from below, 
and odd-order Py approximations converge mono- 
tonically from above, as the order N increases. 

The work on transformations of variational 
principles serves mainly to unify the basic theory. 
However, the bounding property of complementary 
variational principles could have practical application 
to self-adjoint, or extremum, problems. Comple- 
mentary principles could be employed to define upper 
and lower bounds on integral properties or to obtain 
figures of merit for approximate (e.g., synthesis) 
solutions. 


SUMMARY 


Flux synthesis methods provide a powerful means 
for constructing neutron-flux distributions in a degree 
of detail which is impractical, if not impossible, by 
direct solution methods. It is paradoxical that the 
essential feature of flux synthesis methods, the combi- 
nation of detailed lower dimensional solutions to 
obtain a detailed higher dimensional solution, is 
simultaneously their principal advantage and disad- 
vantage. For it can be said fairly of synthesis methods 
that their strongest point is that one can incorporate 
physical insight to reduce an otherwise intractable 
problem to manageable proportions, and, conversely, 
their weakest point is that one must. This, obviously, 
strikes different people in different ways. 

The power and economy of flux synthesis methods 
have been amply demonstrated by numerous applica- 
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tions. to realistic reactor models. These applications 


have 


also. illustrated the need to exercise some 


ingenuity in the selection of trial functions, if accurate 
solutions are required for complex models in which the 
flux is highly nonseparable, and the need to be wary of 
the possibility of anomalous solutions. In the hands of 
an experienced user, flux synthesis methods constitute 
one of the most powerful tools available for the 
analysis of realistic reactor models. 
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dilution capacity, and the discharge structure should be 
designed to promote rapid dilution of the heated effluent by 
the cooler receiving waters. Intake structures should be 
designed to provide for low intake velocities and to provide 
an effective means of diverting fish and other aquatic 
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Test (LOFT) pressurized-water reactor [55 MW(t)] with a 
5.5-ft-long core compares with that of a large four-loop 
pressurized-water reactor [3000 MW(t)] with a 12-ft-long 
core under loss-of-coolant accident conditions, The analysis 
employed the thermal—hydraulic computer codes RELAP3 
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and THETAI-B to determine the hydraulic response of the 
primary systems during blowdown and the thermal responses 
of the cores during blowdown. An additional analysis was 
performed to assess the core thermal response during ECC 
reflooding. On the basis of the results of the analyses, the 
conclusion was reached that a 5.5-ft core operating in the 
LOFT system can produce thermal—hydraulic behavior 
tepresentative of that calculated to occur in a large four-loop 
PWR with a 12-ft core during a double-ended inlet pipe 
break, The inlet pipe break analyzed is generally considered 
to result in the highest core temperatures for large PWRs. 
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Available as ANCR-1053. Date of report, April 1972. 


INDEXED BIBLIOGRAPHY OF THERMAL 


EFFECTS LITERATURE—2 

Authors: J. G. Morgan and C, C, Coutant, Oak Ridge National 
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Abstract: The second volume in the series of indexed bibli- 
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An international symposium sponsored by om 


U.S. Atomic Energy Commission and 


State University of New York at Albany Series 
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EDITORS 


James W. Corbett and Louis C. Ianniello 


884 pages, 6 by 9, paperback 
Library of Congress Catalog Number: 72-600048 


This book summarizes the remarkable progress achieved in just a few years on the important 
problem of voids and swelling in metals subjected to high-temperature irradiation. This swelling 
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Published as No. 26 in the AEC Symposium Series, this book is available as CONF-710601 for $9.00 from 


National Technical Information Service 
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Sodium Technology 1962—1971 (T1ID-3334 Part II) 
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Prepared by the USAEC Technical Information Center for tk: 
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covered is on sodium as applied to reactor-coolant technology. 
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